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Executive  Summary 


Soil  vapor  extraction  (SVE)  is  a  prevalent  remediation  approach  for  volatile  contaminants  in  the  vadose 
zone.  A  diminishing  rate  of  contaminant  extraction  over  time  is  typically  observed  due  to  1)  diminishing 
contaminant  mass,  and/or  2)  slow  rates  of  contaminant  removal  from  low-permeability  zones.  After  a 
SVE  system  begins  to  show  indications  of  diminishing  contaminant  removal  rate,  SVE  performance 
needs  to  be  evaluated  to  determine  whether  the  system  should  be  optimized,  terminated,  or  transitioned  to 
another  technology  to  replace  or  augment  SVE.  For  many  sites,  part  of  this  evaluation  will  include 
consideration  of  vapor  intrusion  into  buildings  and  quantification  of  the  amount  of  source  remediation 
necessary  to  meet  vapor  intrusion  protection  goals.  This  document  presents  a  tool  for  estimating  vapor 
intrusion  into  buildings  that  is  a  result  of  unsaturated  (vadose)  zone  contaminant  sources.  The  tool  builds 
on,  and  is  related  to,  guidance  for  evaluation  of  SVE  performance  relative  to  the  impact  of  a  vadose  zone 
source  on  groundwater  concentrations. 

The  tool  described  here  is  focused  on  vapor-phase  diffusion  from  the  current  vadose  zone  source,  and  is 
not  configured  to  address  vapors  that  arise  either  solely  from  groundwater  plumes  (i.e.,  where 
contamination  was  released  directly  to  groundwater  or  where  a  vadose  zone  source  is  not  currently 
present)  or  from  groundwater  contamination  resulting  from  a  historically  stronger  vadose  zone  source. 
The  spreadsheet  tool  uses  a  quantitative  estimation  approach  that  incorporates  vadose  zone  source 
strength  and  location,  vadose  zone  transport,  and  a  model  for  estimating  movement  of  soil-gas  vapor 
contamination  into  buildings.  The  tool  may  be  appropriate  for  a  broad  range  of  site  conditions,  though 
the  user  will  need  to  assess  the  suitability  of  the  underlying  assumptions  used  in  the  estimation  process 
when  applying  the  tool.  In  some  cases  (such  as  significant  vadose  zone  heterogeneities  or  high 
recharge/moisture  content),  site-specific  modeling  may  be  appropriate  for  improving  estimates. 

The  tool  described  here  incorporates  the  results  of  three-dimensional  multiphase  transport  simulations 
predicting  the  vapor  concentration  distribution  from  a  vadose  zone  source.  A  large  number  of  simulations 
were  conducted  to  develop  a  data  set  that  can  be  queried  based  on  the  specific  configuration  of  the  user’s 
site.  If  the  user’s  site  configuration  does  not  exactly  match  one  of  the  existing  simulation  configurations, 
interpolation  is  used  to  generate  an  appropriate  vapor  concentration  estimate.  With  this  approach,  the 
user  has  access  to  results  from  transport  simulations  without  expending  the  time  and  effort  to  conduct 
simulations  themselves.  This  tool  configuration  also  enables  the  user  to  rapidly  conduct  multiple 
evaluations  using  variations  in  input  parameters,  thereby  gaining  insight  into  the  sensitivity  of  the  results 
to  variations  in  these  parameters. 
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1.0  Introduction 


Soil  vapor  extraction  (SVE)  is  a  prevalent  remediation  approach  for  volatile  contamination  in  the  vadose 
(unsaturated)  zone  of  the  subsurface.  A  diminishing  rate  of  contaminant  extraction  over  time  is  typically 
observed  due  to  1)  diminishing  contaminant  mass  in  the  vadose  zone  source,  and/or  2)  slow  rates  of 
contaminant  removal  from  low-permeability  zones.  After  a  SVE  system  begins  to  show  indications  of 
diminishing  contaminant  removal  rate,  SVE  performance  needs  to  be  evaluated  to  determine  whether  the 
system  should  be  optimized,  terminated,  or  transitioned  to  another  technology  to  replace  or  augment  SVE. 
For  many  sites,  part  of  this  evaluation  will  include  consideration  of  vapor  intrusion  into  buildings  and 
quantification  of  the  amount  of  source  remediation  necessary  to  meet  vapor  intrusion  protection  goals. 

This  document  presents  an  approach  and  a  tool  for  estimating  vapor  intrusion  resulting  from  vadose  zone 
contaminant  sources,  as  depicted  in  Figure  1 .  The  tool  may  be  appropriate  for  readily  obtaining  estimates 
for  a  broad  range  of  site  conditions.  However,  the  user  needs  to  consider  the  underlying  assumptions  and 
limitations  (see  Section  2.5)  of  the  estimation  process  when  applying  the  tool.  The  tool  described  here  is 
focused  on  vapor-phase  diffusion  from  the  current  vadose  zone  source,  and  is  not  configured  to  address 
vapors  that  arise  either  solely  from  groundwater  plumes  ( i.  e. ,  where  contamination  was  released  directly 
to  groundwater  or  where  a  vadose  zone  source  is  not  currently  present)  or  from  groundwater 
contamination  resulting  from  a  historically  stronger  vadose  zone  source.  Further,  the  nature  of  the 
underlying  generalized  conceptual  site  model  necessitates  various  assumptions  and  bounds  for  site 
configuration.  Thus,  a  site-specific  model  may  be  useful  for  detailed  analysis  of  some  sites. 


Commercial/Industrial  Residential 


Groundwater  Plume 


Groundwater 

Figure  1.  Conceptual  representation  of  the  processes  relevant  to  a  constant  vadose  zone  source  at  a 

steady-state  equilibrium.  Contaminant  vapor  diffuses  out  from  the  source  and  may  travel 
within  the  vadose  zone  or  leave  the  vadose  zone  1)  at  the  ground  surface,  2)  via  vapor 
intrusion  into  buildings,  or  3)  by  mass  transfer  into  groundwater  (near  the  source  or  at 
distance  from  the  source).  Groundwater  contaminated  by  the  vadose  zone  source  may  travel 
downgradient  by  advection/dispersion,  and  contamination  could  transfer  back  into  the 
vadose  zone  (depending  on  the  chemical  properties  and  concentration  gradients). 
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When  vapor-phase  transport  is  the  major  component  of  the  overall  contaminant  fate  and  transport  from  a 
vadose  zone  source,  the  contaminant  concentration  expected  at  a  given  subsurface  location  is  controlled 
by  a  limited  set  of  parameters  describing  source  and  subsurface  characteristics.  Thus,  it  is  possible  to  pre¬ 
model  contaminant  transport  for  a  range  of  conditions  and  to  estimate  the  results  for  a  specific  waste  site 
based  on  these  pre-modeled  scenarios.  This  approach  has  previously  been  applied  in  the  Soil  Vapor 
Extraction  Endstate  Tool  (SVEET)  software  [Truex  et  al.,  2013;  Oostrom  et  al.,  2014],  The  tool 
described  herein  incorporates  the  results  of  three-dimensional  multiphase  transport  simulations  (described 
in  detail  by  Truex  et  al.  [2013]  and  Oostrom  et  al.  [2014])  predicting  the  steady-state  vapor  concentration 
distribution  resulting  from  a  vadose  zone  source.  A  large  number  of  simulations  were  conducted  to 
develop  a  data  set  that  can  be  queried  based  on  the  specific  configuration  of  the  user’s  site.  If  the  user’s 
site  configuration  does  not  exactly  match  one  of  the  existing  simulation  configurations,  interpolation  and 
scaling  (e.g.,  to  account  for  differences  in  Henry’s  Law  constants)  are  used  to  generate  an  appropriate 
vapor  concentration  estimate.  With  this  approach,  the  user  has  access  to  results  from  transport 
simulations  without  expending  the  time  and  effort  to  configure  a  model  and  conduct  simulations 
themselves.  This  pre-modeled  scenario  approach  allows  the  user  to  rapidly  conduct  multiple  evaluations 
using  variations  in  input  parameters,  thereby  gaining  insight  into  the  sensitivity  of  the  results  to  variations 
in  these  parameters. 

The  simulations  conducted  for  the  pre-modeled  scenarios  predict  soil  gas  concentrations  for  a  range  of 
subsurface  conditions  where  transport  from  the  vadose  zone  source  occurs  predominantly  by  vapor-phase 
diffusion.  That  is,  the  subsurface  vadose  zone  is  relatively  dry  and  there  is  low  recharge  (infiltration) 
through  the  vadose  zone.  The  pre-modeled  scenario  simulations  included  three-dimensional  contaminant 
movement  from  the  vadose  zone  source  (both  by  diffusion  and  along  with  recharge).  Vapor-phase 
contamination  could  exit  the  vadose  zone  at  the  ground  surface  or  could  partition  into  the  groundwater. 
Contamination  in  the  groundwater,  from  vapor-phase  partitioning  or  recharge  influx,  migrated 
downgradient  with  the  groundwater,  but  without  any  adsoiption  or  bio/geochemical  reactions.  Back- 
transfer  of  contamination  from  groundwater  to  the  vadose  zone  at  downgradient  locations  could  occur, 
depending  on  the  concentration  gradient. 

The  pre-modeled  simulations  represent  neither  historical  groundwater  contamination  nor  groundwater 
contamination  in  the  absence  of  the  current  vadose  zone  source.  That  is,  the  simulations  start  with  clean 
groundwater  and  only  represent  groundwater  contamination  resulting  from  the  defined  vadose  zone 
source.  The  intent  of  this  approach  was  to  allow  evaluation  of  whether  the  residual  source  remaining  after 
termination  of  SVE  is  predicted  to  result  in  a  vapor  intrusion  issue  (assuming  sufficient  time  to  produce  a 
steady-state  vapor  plume).  Such  evaluation  would  support  decisions  regarding  optimization,  termination, 
or  transition  of  the  SVE  remedy.  However,  the  evaluation  approach  described  here  does  not  provide 
insight  into  vapor  intrusion  resulting  from  pre-existing  groundwater  contamination  (i.e.,  that  does  not 
arise  from  the  current  defined  vadose  zone  source). 

The  estimated  soil  gas  contaminant  concentration  obtained  from  the  pre-modeled  scenarios  for  a  building 
location  of  interest  is  combined  with  an  analytical  solution  adapted  from  Johnson  and  Ettinger  [1991]  to 
estimate  the  vapor  intrusion  into  the  building.  This  is  a  different  approach,  in  that  the  soil-gas 
contaminant  concentration  below  the  building  is  determined  independently  and  then  the  entry  rate  of 
contaminant  into  the  building  is  equated  to  the  rate  of  contaminant  removal  via  ventilation.  This  allows 
calculation  of  the  indoor  air  concentration.  The  Johnson  and  Ettinger  [1991]  model  included  equations 
for  diffusive  transport  of  contaminant  in  vadose  zone  soil  gas  from  a  source  (nominally  groundwater 
below  the  building)  and  calculated  an  attenuation  factor  representing  the  ratio  of  indoor  air  concentration 
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to  the  source  concentration.  With  the  independent  calculation  of  concentration  in  the  soil  gas  below  the 
building  applied  in  the  approach  described  here,  there  is  not  a  linkage  between  diffusive  flux  and  rate  of 
vapor  intrusion.  That  is,  the  advective  flow  into  the  building  (driven  by  the  building  under-pressure  and 
air-exchange  with  the  outside)  does  not  influence  the  contaminant  transport  through  the  vadose  zone.  As 
noted  by  the  EPA  [2004],  the  zone  of  influence  for  advective  contaminant  flow  is  expected  to  be  limited 
to  a  region  adjacent  to  the  building  foundation. 

The  overall  vapor  intrusion  estimation  approach  consists  of  four  steps:  1)  defining  site-specific  inputs, 
2)  interpolating  between  pre-modeled  scenario  results  based  on  parameters  that  have  nonlinear  impacts  on 
the  soil  gas  contaminant  concentration,  3)  scaling  the  interpolated  result  based  on  parameters  that  have 
linear  impacts  on  the  soil  gas  contaminant  concentration,  and  4)  estimating  the  contaminant  concentration 
in  a  building  due  to  vapor  intrusion.  The  first  three  steps  result  in  an  estimate  of  the  steady-state  soil-gas 
contaminant  concentration  at  a  depth  and  lateral  distance  representing  the  location  just  beneath  the 
building  of  interest.  The  soil  gas  concentration  is  then  used  as  part  of  the  input  to  the  adapted  Johnson 
and  Ettinger  [1991]  vapor  intrusion  calculations.  For  user  convenience,  the  calculational  steps  are 
incorporated  into  the  Vapor  Intrusion  Estimation  Tool  for  an  Unsaturated-zone  Source  (VIETUS) 
software,  which  is  a  spreadsheet-based  (Microsoft  Excel)  tool. 

This  report  describes  the  calculational  steps  and  basis  for  the  approach,  provides  an  example  to  illustrate 
the  calculations,  and  provides  a  guide  for  using  the  VIETUS  spreadsheet  tool. 
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2.0  Calculational  Steps 


The  approach  for  estimating  a  soil-gas  contaminant  concentration  in  the  vadose  zone  and  the 
corresponding  gas  concentration  in  a  building  due  to  vapor  intrusion  is  based  on  a  generalized  conceptual 
site  model  and,  as  mentioned  previously,  is  comprised  of  several  calculational  steps. 

The  generalized  conceptual  model  describes  a  site  with  a  vadose  zone  source  that  remains  constant  over 
time,  resulting  in  a  steady-state  distribution  of  vapor-phase  contamination  in  the  vadose  zone.  Vapor- 
phase  transport  by  diffusion  from  the  source  zone  is  taken  as  the  dominant  transport  mechanism.  The 
nature  of  the  porous  media  means  that  the  effective  gas  diffusion  in  the  subsurface  is  less  than  that  of 
diffusion  in  air  (because  of  the  pore  structure  of  the  soil  and  a  portion  of  the  pore  space  being  filled  by 
moisture).  Thus,  transport  predominantly  by  diffusion  implies  low  moisture  content  and  low  recharge 
conditions.  Because  the  approach  focuses  on  steady-state  conditions,  moderate  heterogeneities  have 
limited  impact  on  the  vapor-phase  concentration  distribution  [Truex  et  al.,  2013]  and  a  homogenous 
subsurface  approximation  can  be  used  for  the  conceptual  model.  However,  if  a  site  has  significant 
heterogeneities,  site-specific  modeling  may  be  appropriate  for  improving  estimates.  Higher  recharge  and 
wetter  vadose  zone  environments  lead  to  additional  complexity  in  contaminant  transport,  especially  with 
respect  to  interaction  with  the  groundwater.  Such  conditions  are  outside  the  range  of  the  simulations 
conducted  to  support  the  calculational  approach  and  would  also  require  site-specific  modeling. 

The  generalized  conceptual  model  includes  contaminant  partitioning  into  and  out  of  the  groundwater 
beneath  the  source  area  because  this  is  an  important  boundary  condition  for  vapor-phase  transport. 
Because  the  approach  is  intended  to  supply  information  to  support  a  decision  about  SVE  disposition,  the 
calculations  focus  solely  on  the  contribution  to  vapor  intrusion  from  the  current  vadose  zone  source. 
Specifically,  the  only  groundwater  contamination  contributions  included  are  those  directly  resulting  from 
the  current  vadose  zone  source  conditions. 

The  user  should  consider  that  the  generalized  conceptual  model  and  associated  calculation  approach 
include  simplifications  of  actual  site  conditions.  The  approach,  however,  enables  a  user  to  relate 
contaminant  mass  discharge  (or  source  concentration)  from  a  vadose  zone  source  to  a  predicted  steady- 
state  vapor-phase  concentration  at  a  specified  location  for  a  depth  underneath  the  building  of  interest. 
The  soil  gas  concentration  can  then  be  used  in  estimating  potential  vapor  intrusion.  The  user  should 
consider  estimate  uncertainties  in  light  of  the  site  conditions  to  interpret  whether  the  estimated  vapor 
intrusion  concentrations  are  acceptable  or  not. 

The  generalized  conceptual  site  model,  depicted  in  Figure  2,  provides  a  method  for  specifying  the  key 
input  parameters  in  a  framework  that  can  be  adapted  for  a  variety  of  sites.  The  generalized  conceptual 
site  model  represents  the  source  area  as  a  zone  with  specified  dimensions  (i.e.,  areal  footprint,  thickness, 
and  depth  within  the  vadose  zone)  and  contaminant  “strength.”  The  subsurface  is  assumed  to  be  a 
uniform,  homogeneous  material  with  specified  characteristics  (e.g.,  temperature,  moisture  content)  and 
fixed  values  of  porosity  (0.3),  bulk  density  (1.855  g/mL),  and  recharge  (0.4  cm/yr).  Groundwater  beneath 
the  vadose  zone  flows  at  a  specified  rate  (Darcy  velocity),  with  the  orientation  of  the  conceptual  model 
such  that  groundwater  flow  is  parallel  to  the  x-axis  (left  to  right  in  the  figure).  The  model  is  configured  to 
maximize  partitioning  across  the  water  table,  thus  assumes  a  minimal  capillary  fringe.  Contaminants  can 
partition  into  the  groundwater,  move  along  with  the  groundwater  flow,  and,  if  conditions  are  appropriate, 
partition  back  into  the  vadose  zone.  A  single  building  is  assumed  to  be  located  at  a  specified  lateral 
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distance  from  the  center  of  the  source  area,  not  necessarily  on  the  axis  of  the  model.  Characteristics 
pertaining  to  vapor  intrusion  through  the  foundation  (e.g.,  foundation  thickness,  ventilation  rate,  area  of 
cracks,  etc.)  must  be  defined  for  the  building. 


0~C’  Anlet’  Qs>  Qb<  Qc.total’  Pc, bulk'  bJ< 

(Wb,  Lb.  AP.  n,  Wc,  Ac,  xc>-  „ 


Figure  2.  Conceptual  framework  for  estimating  the  impact  of  a  vadose  zone  contaminant  source  on  soil 
gas  concentrations  and  vapor  intrusion  into  a  building;  Tables  1  and  2  describe  the  variables. 


If  the  framework  of  the  generalized  conceptual  site  model,  including  the  valid  ranges  of  parameter  values 
(Section  2.1)  and  the  simplifying  assumptions  (Section  2.5),  is  not  appropriate  for  a  given  site,  then  an 
alternative  calculational  approach  would  be  required.  For  example,  sites  with  significant  hetero¬ 
geneity/layering,  widely  dispersed  contaminant  sources,  significant  contaminant  sources  in  the 
groundwater,  significantly  different  moisture  content,  or  significantly  different  recharge  conditions  may 
need  to  use  a  site-specific  numerical  model  (e.g.,  STOMP  [White  and  Oostrom,  2006])  for  an  improved 
estimate  of  the  soil  gas  concentration  distribution. 

While  each  site  is  unique,  there  are  three  general  categories  of  vadose  zone  contaminant  sources,  as 
depicted  in  Figure  3.  These  categories  reflect  common  types  of  subsurface  heterogeneities  and  are  based 
on  the  premise  that  vadose  zone  contaminant  sources  are  persistent  because  they  reside  within  lower 
permeability  zones  or  areas  that  are  poorly  swept  by  a  SVE  system  (e.g.,  high  moisture  zones).  For  many 
sites,  persistent  vadose  zone  contaminant  sources  would  be  expected  to  exist  within  a  localized  portion  of 
the  vadose  zone  as  a  remnant  of  contaminant  transport  pathways  that  occurred  during  the  waste  disposal 
period.  For  sites  with  widely  dispersed  remaining  sources  and  unique  subsurface  features,  a  site-specific 
analysis  may  be  more  appropriate.  However,  for  sites  where  the  remaining  sources  are  more  localized, 
the  generalized  conceptual  model  may  be  representative  of  the  source  area(s). 
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Figure  3.  Categories  of  conceptual  site  models  for  persistent  vadose  zone  contamination  sources.  A) 
homogenous  subsurface  with  residual  source;  B)  simple  layered  subsurface  with  source 
accumulation  in  low  permeability  material;  and  C)  multiple  layers  or  lenses  in  the  subsurface 
with  source  accumulation  in  low  permeability  material.  Red  =  contamination  source; 
blue  =  groundwater;  tan  =  high  permeability;  brown  =  low  permeability;  black  =  waste 
disposal  site.  The  dashed  lines  show  zones  where,  over  time,  vapor  concentrations  will 
nominally  equilibrate  to  a  composite  source  concentration  as  a  result  of  diffusion  from  the 
specific  source  zones. 


While  the  nature  of  the  subsurface  is  significantly  different  for  the  three  conceptual  model  categories 
shown  in  Figure  3,  the  vadose  zone  source  may  be  approximated  in  each  case  as  a  single  zone  of  specified 
dimension  and  concentration/source  strength  (Figure  2).  This  generalized  source  is  most  apparent  for 
Figure  3(B),  where  remaining  contaminants  reside  primarily  in  a  distinct  low-permeability  layer. 
However,  in  both  Figure  3(A)  and  3(C),  the  contamination  resides  within  a  defined  volume  of  the  vadose 
zone  and  a  composite  concentration  can  be  assigned.  That  is,  for  the  scenarios  shown  in  Figures  3(A)  and 
3(C),  vapors  may  emanate  from  smaller  distinct  portions  of  the  composite  source  zone,  but  in  the  long 
term,  vapor  concentrations  will  equilibrate  in  some  zone  between  and  around  these  distinct  sources  (as 
represented  by  the  dashed  outlines  on  the  figures).  The  result  is  an  effective  composite  source  zone  that 
acts  as  the  source  concentration  that  drives  diffusion  out  into  the  vadose  zone  and  eventually  to  the 
surface  or  the  groundwater.  Aqueous  recharge  moving  through  this  zone  will  become  contaminated. 
However,  the  concentration  of  the  pore  water  at  the  water  table  will  be  in  equilibrium  with  the  vapor 
concentration  at  that  location.  Thus,  for  sites  with  relatively  low  recharge  rates,  the  composite  vapor 
concentration  will  also  be  a  primary  factor  in  the  long-term  aqueous  phase  contaminant  discharge  to  the 
groundwater. 


A  site  may  also  have  heterogeneities  outside  of  the  source  zone  in  the  form  of  distinct  layers  or  lenses 
with  lower-permeability  and/or  higher  moisture  content  than  the  bulk  vadose  zone  materials.  However,  in 
Truex  et  al.  (2013;  in  their  Appendix  E)  conducted  simulations  to  evaluate  long-term  vapor  transport  in 
the  presence  of  lower  permeability  layers  with  moderately  higher  moisture  content  than  the  bulk  vadose 
zone.  Although  the  effective  diffusion  coefficient  is  lower  for  the  increased  moisture  content  in  the  lower 
permeability  layers,  the  subsurface  vapor-phase  contaminant  concentration  distribution  was  only 
minimally  affected  by  such  layers  when  looking  at  long-term  steady-state  conditions.  Thus,  for  moderate 
subsurface  heterogeneities,  the  use  of  the  homogenous  approximation  in  the  generalized  conceptual 
model  may  be  appropriate.  Users  will  need  to  consider  the  level  of  heterogeneity  at  their  site  and  assess 
the  appropriateness  of  the  generalized  conceptual  model,  using  site-specific  modeling  for  the  vapor 
intrusion  assessment  if  the  generalized  conceptual  model  is  insufficient  for  their  site. 
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The  four  steps  of  the  calculational  approach  implemented  in  the  VIETUS  software  are  shown  in  Figure  4, 
and  are  discussed  in  detail  in  the  sections  below.  Section  3.0  walks  through  an  example  calculation  to 
further  illustrate  these  steps. 


Figure  4.  Flow  chart  of  the  steps  involved  in  the  process  for  estimating  soil  gas  contaminant 

concentration  and  building  contaminant  concentration  due  to  vapor  intrusion;  Sections  2. 1 
through  0  correspond  to  these  four  steps. 


2.1  Compile  Inputs  for  the  Conceptual  Site  Model  Framework 

The  framework  of  the  generalized  conceptual  site  model  (Figure  2)  requires  a  set  of  inputs  to  describe 
each  scenario  of  interest.  The  input  parameters  listed  in  Table  1  describe  the  source-  and  transport-related 
characteristics  that  must  be  compiled  (with  vapor  intrusion  parameters  being  discussed  later  in  this 
section).  These  inputs  are  used  to  calculate  several  intermediate  parameters.  The  inputs  (or  the 
appropriate  intermediate  parameter)  define  the  interpolation  between  pre-modeled  scenarios  (for 
parameters  exhibiting  a  nonlinear  response  in  the  results)  and  the  linear  scaling  from  the  base  case 
scenario  to  the  site-specific  scenario.  To  avoid  extrapolation  outside  of  the  pre-modeled  scenarios,  the 
input  parameters  are  restricted  to  be  within  the  permissible  ranges  noted  in  Table  1.  Ranges  have  also 
been  defined  for  parameters  that  exhibit  a  linear  or  inverse  linear  relationship  to  contaminant 
concentrations  (including  the  “source  strength”  parameters).  The  source  strength  input  parameters  (Cgs 
and  MSrc)  have  upper  bounds  nominally  based  on  levels  where  density-driven  advection  would  not  affect 
contaminant  migration.  The  range  for  Henry’s  Faw  constant  is  constrained  by  the  choice  of  the 
contaminant  and  the  permissible  range  for  the  subsurface  temperature. 


Table  1.  Parameters  defining  source  and  transport  site  characteristics  used  in  the  estimation  of  soil  gas 
concentration  (shaded  rows  are  calculated  parameters;  rows  without  shading  are  user  inputs) 


Parameter  Name 

Units 

Variable 

Permissible  Range 

Key  Values a 

INPUT  PARAMETERS  USED  to  calculate  other  parameters 

Contaminant 

— 

— 

Carbon  Tetrachloride 
Chloroform 
Dichloromethane 
Chloromethane 
Tetrachloroethene 
Trichloroethene 
Dichloroethene  isomers 
Vinyl  Chloride 

1 .1 .1 .2- Tetrachloroethane 

1 .1 .2.2- Tetrachloroethane 

1.1.1- Trichloroethane 

1 .1 .2- Trichloroethane 

1.1- Dichloroethane 

1 .2- Dichloroethane 
Chloroethane 

Carbon 

Tetrachloride 

Average  subsurface  temperature 

°C 

T 

10-30 

20 

Depth  to  the  top  of  the  source  zone 

m 

LI 

0.5-49 

— 

Vertical  thickness  of  the  source  zone 

m 

z 

1  -30 

— 

Width  for  plan  view  extent  of  source 
(source  is  assumed  square,  thus  width  =  length) 

m 

w 

10-50 

— 

Parameters  having  a  nonlinear  impact  on  concentrations 

Average  gravimetric  moisture  content  in  vadose  zone b 

(mass  of  water  per  mass  of  dry  porous  medium) 

wt% 

UJ 

1  -9 

1, 5,  9 

Vadose  zone  thickness 

m 

VZT 

10-60 

10,  30,  60 

Source  thickness  ratio  (fraction  of  vadose  zone 
thickness  that  is  source  material) 

(STR  =  z/  VZT) 

— 

STR 

o 

1 

o 

cn 

0.1, 0.25,  0.5 

Relative  vertical  source  position  in  the  vadose  zone 
(i.e.,  near  surface,  aquifer  or  middle) 

(RSP  =  LI  / L2  =  LI  / [VZT -  LI -z]) 

— 

RSP 

0.1  -  10 

0.1,  1,  10 

Source  area  as  a  plan-view  areal  extent 

(SA  =  w 2) 

m2 

SA 

100-2500 

100,  400, 
900, 2500 

Groundwater  Darcy  velocity 
(q  is  input  directly,  but  can  be  calculated  from  q  =  Kh-i,  with 

Kh  =  horizontal  hydraulic  conductivity,  i  =  hydraulic  gradient) 

m/d 

q 

0.005-0.3 

0.005,  0.03, 
0.3 

Parameters  having  a  linear  (or  inverse  linear)  relationship  with  soil  gas  concentration 

Henry’s  Law  constant0 

— 

H 

contaminant-specific 

0.89 

Gas-phase  contaminant  concentration  in  the  source01 

ppmv 

Cgs 

0.1  -2000 

159 

Mass  discharge  of  vapor-phase  contamination  from 
the  source01 

g/d 

Msrc 

0.1  -5000 

varies 

a  The  key  values  indicate  either  the  values  used  in  the  pre-modeled  scenario  simulations  (for  parameters  having  nonlinear  relationship)  or 
the  “base  case”  scenario  values  (for  parameters  having  a  linear/inverse  linear  or  source  strength  relationship).  Base  case  recharge  was 
0.4  cm/yr,  but  recharge  is  not  a  site-specific  input.  The  pre-modeled  scenario  simulations  are  discussed  in  more  detail  in  Section  2.2. 
b  The  pre-modeled  scenario  simulations  (Section  2.2)  actually  use  the  residual  water  saturation  (Sr),  which  is  calculated  from  moisture 
content  as:  Sr  =  (cjJ-pbUik)  /  (100-0, otai'Pwater),  where  pbuik  is  the  dry  bulk  density  of  the  porous  media  (1 .855  g/mL),  0,o,ai  is  the  total  porosity 
(0.3),  and  pwater  is  the  (temperature-dependent)  density  of  water.  The  range  for  S,  was  0.05  to  0.55  in  the  simulations,  so  w  is 
constrained  accordingly,  as  discussed  in  the  text  in  below. 

0  The  dimensionless  Henry’s  Law  Constant  is  based  on  the  site-specific  subsurface  temperature  and  contaminant-specific,  temperature- 
dependant  property  correlations  for  vapor  pressure  [Yaws  et  al.,  2009]  and  solubility  [Mackay  et  al. ,  2006:  Yaws,  2003]. 
d  The  user  may  provide  either  CgS  or  MDS,  but  not  both.  The  mass  discharge  “key  values”  are  interpolated  from  the  pre-modeled  scenario 
simulation  results. 
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Note  that,  while  gravimetric  soil  moisture  content  (co)  is  requested  as  an  input  parameter  (Table  1)  for 
user  convenience,  residual  saturation  (Sr)  is  the  actual  parameter  specified  in  the  pre-modeled  scenario 
simulations  (Section  2.2).  Sr  is  calculated  from  co  as  discussed  in  the  footnotes  to  Table  1.  The  key 
values  for  Sr  were  0.05,  0.3,  and  0.55,  which  correspond  to  moisture  content  values  of  0.8071721,  4.843, 
and  8.87889,  respectively,  for  a  subsurface  temperature  of  20  °C.  To  keep  the  information  concise,  the 
permissible  moisture  content  range  is  specified  in  Table  1  and  the  VIETUS  software  as  1  to  9  wt%. 
Values  at  or  above  8.87889  wt%  are  treated  as  Sr  values  of  0.55.  If  the  “Allow  ©  down  to  Sr  =  0.05?” 
option  (discussed  in  Section  4.4)  is  set  to  TRUE,  then  co  values  down  to  0.8071721  wt%  (at  20  °C)  are 
allowed.  Otherwise,  co  values  below  1  wt%  result  in  an  incomplete  calculation  (because  the  value  is 
outside  the  permissible  range).  The  specific  weight  percent  numbers  mentioned  above  vary  slightly  as  a 
function  of  temperature  because  the  density  of  water  is  a  function  of  temperature.  Figure  5  depicts  the 
permissible  range  for  co  and  the  corresponding  values  of  S, .  Note  that  the  numbers  discussed  above 
include  sufficient  significant  digits  to  obtain  precision  in  the  appropriate  conversion  calculations. 


OJ 

Sr 


Only  permitted  when 
the  "Allow  to  down  to 
Sr  =  0.05?"  option 
is  set  to  TRUE 


A 


& 


& 


& 


0.55 


(wt%) 

(-) 


Calculated  values  for  a  temperature  of  2CTC  are 
shown  with  sufficient  significant  digits  to  obtain 
precision  in  the  appropriate  conversion  calculations. 


Figure  5.  Depiction  of  the  permissible  range  for  the  gravimetric  soil  moisture  content  input  parameter 

(co),  the  corresponding  values  of  residual  saturation  (5,),  and  how  the  ranges  are  reconciled  at 
the  lower  and  upper  limits 


Appendix  A  of  Truex  et  al.  [2013]  discusses  approaches  for  determining  or  estimating  certain  information 
about  the  source  area,  including  the  lateral  extent,  location  within  the  vadose  zone,  and  “source  strength.” 
The  source  strength  calculational  input  may  be  represented  as  either  a  constant  vapor-phase  contaminant 
concentration  (Cgs)  or  mass  discharge  of  contaminant  (Msrc),  but  not  both. 

Input  parameters  for  calculating  vapor  intrusion  (migration  of  contaminant  from  soil  gas  into  a  building) 
are  based  on  the  Johnson  and  Ettinger  [1991]  model,  with  numerical  simulation  results  used  for  the  trans¬ 
port  component.  Unlike  the  Johnson  and  Ettinger  [1991]  model,  the  VIETUS  calculational  approach 
estimates  a  vapor-phase  contaminant  concentration  in  the  subsurface  below  the  building,  based  on 
numerical  simulations  for  three-dimensional,  multi-phase  transport  from  the  defined  vadose  zone  source 
area.  The  vadose  zone  source  strength  is  held  constant  and  a  steady-state  vapor-phase  contaminant 
distribution  is  determined  in  the  numerical  simulations.  This  approach  thus  ties  the  independently 
calculated  subsurface  vapor-phase  contaminant  concentration  to  the  contaminant  entry  into  a  building, 
which  is  determined  from  the  Johnson  and  Ettinger  [1991]  model.  This  coupling  of  soil-gas 
concentration  distribution  from  pre-modeled  scenarios  to  the  contaminant  influx  to  the  building  in  the 
VIETUS  software  allows  the  user  to  quickly  set  up  and  assess  multiple  building  scenarios.  However,  the 
results  can  differ  from  the  Johnson  and  Ettinger  [1991]  model.  One  difference  is  that  the  Johnson  and 
Ettinger  [1991]  model  allows  the  vapor-phase  concentration  below  the  building  to  vary  to  achieve 
equilibrium  between  the  diffusive  flux  through  the  soil  and  the  advective/diffusive  flux  into  the  building. 
The  pre-modeled  simulations  do  not  include  advective  movement  of  soil  gas,  either  by  barometric 
pumping  or  driven  by  building  under-pressure.  However,  EPA  [2004]  describes  the  influence  of 
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convection  resulting  from  inflow  to  a  building  as  only  affecting  the  zone  immediately  adjacent  to  the 
foundation.  Still,  VIETUS  results  may  not  fully  account  for  relevant  advection  when  the  building  is 
laterally  situated  near  the  source  area. 

The  vapor  intrusion  parameters  required  for  use  of  the  VIETUS  tool  describe  the  foundation,  gas  flow 
pathways  through  the  foundation,  and  the  building.  Table  2  lists  the  input  parameters  relevant  to  vapor 
intrusion  calculations  in  three  categories.  The  first  group  includes  key  parameters  required  for  the 
calculations,  the  second  group  are  key  parameters  that  can  either  be  directly  specified  or  can  be  calculated 
from  secondary  parameters  (based  on  certain  assumptions),  and  the  third  group  consists  of  those 
secondary  parameters. 


Table  2.  Input  parameters  defining  vapor  intrusion  characteristics 


Parameter  Name 

Units 

Variable 

Permissible  Range 

Key  Vapor  Intrusion  Parameters 

(directly  specified) 

Longitudinal  distance  to  the  building 

(from  the  source  center  to  the  center  of  the  building) 

m 

dx 

-950  -  950  a 

Longitudinal  distance  to  the  building 

(from  the  source  center  to  the  center  of  the  building) 

m 

dy 

0-470 

Depth  of  basement/foundation 

m 

dz 

0.5,  1,  1.5,  2,  2.5,  3,  3.5,  or  4 

Foundation/wall  thickness 

m 

Lc 

Building  ventilation  rate 

m3/d 

Qb 

Key  Vapor  Intrusion  Parameters 

(calculated  from  secondary  parameters  or  directly  specified) 

Gas  flow  rate  into  the  building  from  the  soil 

m3/d 

Qs 

0<Qs<(xc/40)-15b’c 

Area  of  cracks  or  foundation 

m2 

Ajnlet 

varies  w 

Total  porosity  within  the  cracks/foundation 

— 

9c,  total 

0.01  -0.9 

Gas-filled  porosity  within  the  cracks/foundation 

— 

9c,  g 

-  9c,  total 

Secondary  Vapor  Intrusion  Parameters 

Vapor  intrusion  path  type 

— 

— 

Cracks,  Porous  foundation 

Total  area  of  cracks 

m2 

Ac 

<  4  -Xcdz 

Building  under-pressure 

Pa 

AP 

0-20 

Soil  gas  permeability 

m2 

kv 

IE-10  -  IE-16 

Crack  ratio  (area  of  cracks  /  area  of  foundation) 

— 

n 

0.0001  -0.002 

Crack  width 

m 

Wc 

0.0005-0.01 

Building  width 

m 

wb 

>  2  and  Wb-Lb  <  900 

Building  length 

m 

Lb 

>  2  and  Wb-Lb  <  900 

Cumulative  lateral  crack  length 

m 

Xc 

>  Ac  1(4- dz) 

a  Negative  distances  are  located  upgradient  (with  respect  to  groundwater  flow  direction)  from  the  source  area. 
b  The  combination  of  the  foundation/wall  thickness,  Lc,  the  crack/foundation  area,  Ainiet,  the  gas  diffusion  coefficient  for  the 

cracks/foundation,  Dc,  and  the  Qs  value  must  not  cause  the  exponential  function  to  result  in  an  error  in  the  following  expression: 
exp[(Lc-Qs)/(Ainiet-Dc)].  If  the  site  configuration  causes  an  error  in  the  exp()  term,  refer  to  equation  22  of  Johnson  &  Ettinger  [1991], 
0  The  upper  limit  shown  for  Qb  is  a  soft  limit,  with  a  hard  limit  applied  at  a  value  of  80,000  m3/d.  The  upper  limit  shown  for  Qs  is  a 
soft  limit,  with  a  hard  limit  applied  at  a  value  of  0.06-Qb  (which  is  twice  the  observed  top  end  sub-slab  to  indoor  air  attenuation 
factor  [EPA,  2012]). 

d  Clearly,  the  area  of  cracks  cannot  exceed  the  area  of  the  subsurface  foundations.  Typically,  the  area  of  cracks  is  less  than  0.2% 
of  the  area  of  the  subsurface  foundations. 
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The  key  parameters  for  vapor  intrusion  calculations  define  the  location  of  the  building  with  respect  to  the 
vadose  zone  contamination  source  and  gas  flow  into  the  building  (whether  through  a  slab  foundation  or 
into  a  sub-grade  basement).  The  lateral  location  of  the  building  and  depth  of  the  foundation  are  used  to 
determine  the  soil  gas  concentration  underneath  the  building  (based  on  numerical  transport  simulation 
results).  Whether  the  lateral  location  represents  the  center  of  the  building  footprint  or  another  location  of 
interest,  a  single  soil  gas  value  is  used  to  calculate  vapor  intrusion. 

Johnson  and  Ettinger  [1991]  describe  several  scenarios  and  apply  some  simplifying  assumptions  to  arrive 
at  estimates  of  migration  of  contaminant  through  a  foundation.  The  first  scenario  they  considered 
involves  cracks  in  the  foundation,  through  which  gas  (and  contaminant)  flows  because  of  building  under¬ 
pressure.  In  the  crack  scenario,  the  area  of  the  cracks  is  specified  and  cracks  are  assumed  to  be  filled  with 
porous  media  having  the  same  porosity,  bulk  density,  and  moisture  content  as  the  soil  underneath  the 
foundations.  In  addition,  Johnson  and  Ettinger  [1991]  estimate  the  volumetric  flow  rate  of  soil  gas  into 
the  building  based  on  an  idealized  model  of  cracks  being  located  at  floor/wall  seams.  The  other  scenario 
they  considered  is  vapor  intrusion  through  a  porous  foundation,  in  which  case  soil  gas  flows  through 
foundation  pores  across  the  entire  sub-grade  portion  of  the  foundation.  These  scenarios  and  assumptions 
provide  context  for  what  input  parameters  are  needed  for  the  calculations  to  estimate  vapor  intrusion. 
Section  0  discusses  the  vapor  intrusion  parameters  in  more  detail. 

Because  of  the  uncertainty  associated  with  input  parameters,  users  are  encouraged  to  conduct  the 
estimation  process  using  a  sensitivity  analysis  approach.  In  such  an  approach,  multiple  estimates  are 
made  using  appropriate  ranges  for  input  parameter  values  and  the  effect  of  these  variations  on  the 
estimated  soil-gas  contaminant  concentration  can  be  assessed.  The  spreadsheet  tool  (Section  4.0)  allows 
multiple  input  scenarios  to  be  evaluated  as  a  group. 

2.2  Interpolation  from  Pre-Modeled  Scenario  Results  for  Nonlinear 
Variables 

A  total  of  972  pre-modeled  scenarios  were  simulated  with  the  Subsurface  Transport  Over  Multiple  Phases 
(STOMP)  code  [White  and  Oostrom,  2006]  to  assess  the  impact  on  results  when  parameters  that  exhibit  a 
nonlinear  response  are  varied.  Simulations  were  run  for  each  scenario  to  obtain  steady-state 
concentration  distributions  in  the  gas  and  aqueous  phases  throughout  the  computational  domain,  assuming 
a  constant-strength  vadose  zone  source.  These  pre-modeled  scenarios  all  used  “base  case”  values  for 
parameters  exhibiting  a  linear  response,  including  values  of  0.4  cm/yr  of  recharge,  a  carbon  tetrachloride 
source  at  a  soil-gas  concentration  of  1  mg/Lgas  (159  ppmv),  and  a  Henry’s  Law  constant  of  0.89.  The  base 
case  also  used  a  porosity  of  0.3,  a  dry  bulk  density  of  1855  kg/m3,  a  gas  diffusion  coefficient  of  0.0845 
cm2/s,  and  a  liquid  diffusion  coefficient  of  8.62E-6  cm2/s.  Table  1  lists  the  key  values  used  in  the 
simulations  for  each  of  the  parameters  with  a  nonlinear  relationship  to  contaminant  concentrations.  It  is 
the  combinations  of  these  key  parameter  values  that  comprise  the  suite  of  972  pre-modeled  scenarios. 
Key  parameter  values  were  selected  to  allow  linear  interpolation  between  simulation  results  (i.e.,  for 
bounding  key  values)  to  provide  a  reasonable  estimate  of  the  soil-gas  concentration  at  any  given 
parameter  value  (within  the  bounds  of  the  key  values).  Note  that  the  pre-modeled  simulations  used  for 
VIETUS  are  the  same  set  of  simulations  as  used  for  the  SVEET  software  [Truex  et  al.,  2013]. 

The  interpolation  step  in  the  estimation  of  the  soil-gas  contaminant  concentration  is  based  on  a  sequence 
of  data  lookups  and  linear  interpolations  to  find  the  unsealed  gas  concentration  (Cgu)  for  the  site-specific 
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input  parameters.  Interpolation  is  needed  for  parameters  and  when  site-specific  parameter  input  values 
differ  from  key  values.  Interpolated  values  are  calculated  using  Equation  1,  where  P  denotes  the 
parameter  value,  C  is  the  simulation  concentration,  and  the  subscripts  upper  and  lower  represent  the 
known  values  corresponding  to  key  values  above  and  below  the  interpolation  point  of  interest  ( interp ). 


C 


interp 


P  -  p 

interp  lower 

P  -  P 

upper  lower 


(c  ~c  ur 

upper  ^ lower  '  ^ lower 


(1) 


Spatial  interpolation  is  generally  required  between  soil-gas  concentration  results  at  grid-cell  centers  in  the 
numerical  model  to  obtain  a  concentration  result  at  the  site-specific  vertical  and  lateral  location.  Soil  gas 
concentrations  at  the  center  of  grid  cells  were  tabulated  for  each  pre-modeled  scenario  to  allow  spatial 
interpolation  across  the  model  domain  down  to  a  depth  of  four  meters  below  ground  surface.  Thus, 
interpolation  can  be  performed,  as  needed,  to  determine  the  soil  gas  concentration  at  the  vertical  location, 
dz,  and  lateral  location  (specified  by  dx  and  dy)  of  interest. 

A  set  of  26  pre-modeled  scenarios  provide  simulation  soil-gas  concentration  values  ( Cv„„)  for  a  given 
location  (at  dx,  dy,  and  dz)  that  can  be  used  for  parameter-based  interpolation  related  to  the  six  parameters 
that  exhibit  a  nonlinear  contaminant  concentration  response  ( RSP ,  q,  SA,  VZT ,  STR,  and  .S', ;  described  in 
Table  1).  These  64  values  represent  32  pairs  of  concentrations  for  scenarios  with  parameter  key  values 
that  are  both  above  and  below  (or  equal  to)  the  site-specific  parameter  value.  Interpolation  between 
parameter  key  values  proceeds,  first  for  RSP,  then  q,  and  so  on  for  SA,  VZT,  STR,  and  Sr,  to  arrive  at  a 
value  for  Cgu  .  Note  that  .S',  is  the  water  saturation  value  converted  from  the  co  moisture  content  value  (as 
discussed  in  the  footnotes  to  Table  1). 


If  the  source  strength  input  parameter  was  provided  as  a  Msrc  value  (i.e.,  not  a  Cgs  value),  then  a  second 
sequence  of  lookups/interpolations  (in  the  same  order  of  RSP,  q,  SA,  VZT,  STR,  and  Sr)  is  performed  to 
determine  the  simulated  contaminant  mass  discharge  (Ms;m)  corresponding  to  the  input  site  parameters. 
This  MSim  value  is  needed  as  a  linear  scaling  factor.  The  process  for  obtaining  the  interpolated  MSm  mass 
discharge  value  is  the  same  as  for  Cgu ,  except  that  spatial  interpolation  does  not  apply. 

For  convenience,  the  key  VZT  and  STR  parameter  values  are  paired  and  assigned  a  Source  Assessment 
Category  designation,  as  shown  in  Table  3.  It  is  infeasible  to  assemble  the  C«m  data  for  soil  gas 
concentrations  into  discreet  tables  organized  by  these  categories  because  the  soil  gas  data  is  comprised  of 
multiple  full  layers  of  data  and  is  about  75  MB  when  stored  as  an  uncompressed,  plain  text  file. 
However,  the  category  terminology  is  useful  nomenclature,  consistent  with  SVEET  [Truex  et  al.,  2013]. 


Table  3.  Categorization  of  pairs  of  VZT  and  STR  values  to  assist  in  data  organization 


VZT  (m) 

STR  (-) 

10 

30 

60 

0.1 

Category  1 

Category  2 

Category  3 

0.25 

Category  4 

Category  5 

Category  6 

0.5 

Category  7 

Category  8 

Category  9 
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2.3  Scaling  for  Linear  Variables 

The  third  step  in  the  estimation  procedure  is  to  scale  the  Cgu  value  obtained  from  the  interpolation  step  to 
account  for  the  parameters  where  the  soil  gas  contaminant  concentration  varies  linearly  or  inverse  linearly 
with  the  parameter  value.  The  base  case  (key)  values  discussed  above  and  listed  in  Table  1  form  the  basis 
for  the  scaling. 

The  Henry’ s  Law  constant  for  the  site  conditions  (contaminant  and  temperature)  is  required  as  part  of  the 
scaling  process.  In  Table  1  it  was  noted  that  the  Henry’s  Law  constant  can  be  calculated  based  on  the 
site-specific  subsurface  temperature  and  contaminant-specific,  temperature -dependent  property  correla¬ 
tions.  The  Henry’s  Law  constant  and  its  temperature  dependence  have  been  examined  in  a  wide  range  of 
literature  for  contaminants  of  environmental  interest  [e.g.,  Staudinger  and  Roberts,  2001;  Warneck,  2007; 
Chen  et  al.,  2012].  Brennan  et  al.  [1998]  suggest  calculating  the  Henry’s  Law  constant  as  the  ratio  of  the 
vapor  pressure  to  the  water  solubility  as  the  preferred  approach  for  dilute  contaminant  concentrations 
(<0.02  mol  fraction).  Thus,  a  temperature-dependent  Henry’s  Law  constant  can  be  found  using 
temperature -dependent  vapor  pressure  and  water  solubility  values.  However,  the  accuracy  of  this 
approach  depends  on  the  accuracy  of  the  vapor  pressure  and  water  solubility  information. 


The  temperature-dependent  vapor  pressure  correlation  selected  for  use  in  this  work  is  the  Antoine 
correlation  given  in  Equation  2  where  T  is  temperature  in  °C,  Pvap  is  the  vapor  pressure  in  mm  Hg,  and  A, 
B ,  and  C  are  contaminant-specific  correlation  coefficients  [Yaws  et  al.,  2009],  listed  in  Table  4. 


Lo8w(Pvap)  =  A 


B 

T  +  C 


Tis  in  °C  and  P vap  is  in  mm  Hg 


(2) 


Table  4.  Tabulated  correlation  coefficients  for  calculating  vapor  pressure  and  solubility  of  contaminants 
of  interest,  taken  from  Yaws  et  al.  [2009],  Mackay  et  al.  [2006],  and  Yaws  [2003] 


Correl.  Coef.  for  Calculating  Vapor  Pressure 

Correlation  Coefficients  for  Calculating  Solubility 

Contaminant 

Abbrev. 

Contaminant 

Weight 

(g/mol) 

A 

B 

c 

A 

B 

C 

D 

E 

F 

CT 

Carbon  Tetrachloride 

153.823 

7.01144 

1278.54 

232.888 

9.7842  E-2 

-1.4942  E-3 

3.5854  E-5 

2.2775  E-7 

0 

0 

CF 

Chloroform 

119.378 

7.11148 

1232.79 

230.213 

9.9500  E-1 

-1.0531  E-2 

7.9819  E-5 

6.6431  E-7 

0 

0 

DCM 

Dichloromethane 

84.933 

7.11464 

1152.41 

232.442 

1.9610 

-4.4883  E-2 

8.6617  E-4 

4.9463  E-6 

0 

0 

CM 

Chloromethane 

50.488 

6.99771 

870.17 

235.586 

1.4019 

-6.3562  E-2 

1.7198  E-3 

-2.8262  E-5 

2.5268  E-7 

-9.3470  E-10 

PCE 

Tetrachloroethene 

165.833 

7.06892 

1458.45 

226.986 

2.6479  E-1 

-1.5487  E-3 

2.4477  E-6 

0 

0 

0 

TCE 

Trichloroethene 

131.388 

6.87981 

1157.83 

202.58 

1 .4049 

-8.2223  E-3 

1.3218  E-5 

0 

0 

0 

1,1  -DCE 

1,1-Dichloroethene 

96.943 

7.21678 

1181.12 

240.84 

6.2741 

-3.8257  E-2 

6.0461  E-5 

0 

0 

0 

cDCE 

cis-1 ,2-Dichloroethene 

96.943 

7.21953 

1290.28 

236.887 

2.7735  E+1 

-1.7832  E-1 

2.9328  E-4 

0 

0 

0 

tDCE 

trans-1 ,2-Dichloroethene 

96.943 

7.21356 

1244.35 

239.497 

7.8039 

-4.5457  E-2 

6.9676  E-4 

0 

0 

0 

VC 

Vinyl  Chloride 

62.498 

6.91423 

911.15 

239.8 

-4.6670  E+2 

1 .3479  E+4 

1 .8523  E+2 

-1.3236  E-1 

0 

0 

1,1,1,2-TeCA 

1,1,1 ,2-Tetrachloroethane 

167.849 

7.03897 

1467.16 

222.34 

2.1790 

-1.3966  E-3 

2.9328  E-5 

0 

0 

0 

1,1,2,2-TeCA 

1 , 1 ,2,2-Tetrachloroethane 

167.849 

6.91043 

1378.88 

197.086 

4.8798 

-3.0937  E-3 

5.2051  E-5 

0 

0 

0 

1,1,1-TCA 

1,1,1-Trichloroethane 

133.404 

7.00718 

1253.2 

229.624 

1.0909 

-6.5278  E-3 

1.1075  E-5 

0 

0 

0 

1,1,2-TCA 

1,1,2-Trichloroethane 

133.404 

7.14357 

1457.65 

228.099 

2.8980 

-1.8585  E-2 

3.4896  E-5 

0 

0 

0 

1,1 -DCA 

1,1-Dichloroethane 

98.959 

7.18316 

1269.43 

237.755 

9.4136 

-5.7249  E-2 

9.1784  E-5 

0 

0 

0 

1 ,2- DCA 

1,2-Dichloroethane 

98.959 

7.29525 

1407.85 

235.48 

1.7915  E+1 

-1.1684  E-1 

2.0003  E-4 

0 

0 

0 

CE 

Chloroethane 

64.514 

7.13047 

1097.6 

246.009 

5.8420  E-1 

-1.6863  E-3 

9.3949  E-5 

-2.5316  E-6 

0 

0 

For  all  potential  contaminants  of  interest  except  vinyl  chloride,  a  polynomial  correlation  [Mackay  et  al., 
2006]  is  used  to  obtain  the  temperature -dependent  solubility,  as  shown  in  Equation  3.  Here,  xp  is  the  mass 
fraction  in  weight  percent,  T  is  temperature  in  K  (unless  the  contaminant  is  CT,  CF,  DCM,  CM,  or  CE,  in 
which  case  the  temperature  is  in  °C),  and  A,  B,  C,  D,  E,  and  F  are  tabulated  contaminant-specific 
correlation  coefficients  [Mackay  et  al.,  2006],  listed  in  Table  4.  Where  fewer  than  six  coefficients  are 
given  by  Mackay  et  al.  [2006],  the  missing  coefficients  are  set  to  0.0.  Equation  4  [Yaws,  2003]  is  used  to 
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calculate  the  solubility  for  VC,  where  x/  is  the  mass  fraction,  T  is  the  temperature  in  K,  and  A,  B,  C,  and  I) 
are  tabulated  contaminant-specific  correlation  coefficients  [Yaws,  2003]  (Table  4).  Mass  fraction  values 
are  converted  to  mole  fraction,  x,  in  Equation  5  by  multiplying  by  the  ratio  of  the  molecular  weight  of 
water  ( MWw ,  18.01528  g/mol)  to  the  molecular  weight  of  the  contaminant  (MW/,  g/mol),  also  in  Table  4. 

xp  =  A  +  B  T  +  C-T2  +  D-T3  +  E-T4  +  F-T5  T  is  in  K,  except  use  °C  for  CT,  CF,  DCM,  CM,  or  CE  (  3  ) 

Logio(Xf)  =  A  +  B/T  +  C-  Logio(T)  +  D-  T  T  is  in  K  (  4 ) 

MW  MWw  (5) 

100  MW,  f  MW , 


As  discussed  above,  the  dimensionless  Henry’s  Law  constant  itself  is  calculated  from  the  ratio  of  the 
vapor  pressure  ( Pvap )  to  the  mole  fraction  (x),  with  appropriate  conversions  from  units  of  atm/mol  fraction 
to  units  of  concentration  per  concentration  (i.e.,  dimensionless).  See  Sander  [1999]  for  a  thorough 
exposition  on  distinguishing  between  multiple  representations  of  the  units  for  the  Henry’s  Law  constant. 
Equation  6  shows  the  calculation  for  the  unitless  Henry’s  Law  constant  (//),  where  pw  is  the  density  of 
water  (g/mL),  Rgas  is  the  gas  constant  (0.08205746  L-atm-  K  '■  mol1),  T  is  the  average  subsurface 
temperature  (K),  and  other  quantities  were  defined  above. 


P 

_  yap 


MW „ 


1  L 


1  atm 


Pw  '  Rgas '  T  1,1  R  760  mmHg 


T  is  in  K 


(6) 


The  unsealed  soil  gas  concentration,  Cgu,  can  be  scaled  to  a  final  concentration  based  on  the  site-specific 
values  of  the  Henry’s  Law  constant  ( H )  and  the  source  strength.  If  the  input  included  a  value  for  Cgs, 
then  Equation  7  is  used  to  calculate  the  final  estimated  soil  gas  contaminant  concentration,  Cg,  at  the 
location  directly  below  the  building  foundation.  If  the  input  included  Msrc ,  then  Equation  8  is  used  to 
calculate  Cg . 


C 


=  c  °-89 

su  H  159.0 


c=c 


0.89  M. 


gu 


H  Msim 


(7) 

(8) 


Note  that  scaling  to  account  for  differences  in  contaminant  gas  diffusion  coefficients  is  not  performed 
because  the  differences  in  coefficient  values  are  small.  For  the  contaminant  species  in  VIETUS,  the  gas 
diffusion  coefficients  average  0.09  ±  0.015  cm2/s  and  the  range  of  the  differences  in  gas  diffusion 
coefficients  relative  to  the  base  case  varies  from  -18%  to  +48%. 


2.4  Estimation  of  Vapor  intrusion 

Vapor  intrusion  is  calculated  based  on  the  approach  described  by  Johnson  and  Ettinger  [1991],  with  the 
adaptation  that  the  soil  gas  concentration  underneath  the  building  foundation,  Cg ,  is  determined  from  the 
interpolation  and  scaling  steps  described  above.  Johnson  and  Ettinger  specify  equations  for  the  entry  rate 
of  vapor  contamination  into  a  building  and  removal  of  a  contaminant  as  a  result  of  building  ventilation  (in 
their  equations  14  and  17,  respectively).  Rearranging  their  equations  to  solve  for  the  concentration  in  the 
building,  Cb,  results  in  Equation  9  below.  Several  key  parameters  (Table  2)  and  one  chemical  property 
are  required  as  inputs  to  this  equation. 
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Estimates  of  appropriate  values  for  the  Equation  9  input  parameters  (Qb,  Ai„iet,  Lc,  and  Qs)  can  be 
obtained  from  relevant  calculations  (some  of  which  are  left  for  the  user  to  complete  external  to  the 
VIETUS  software  in  the  interest  of  streamlined  inputs)  and  from  guidance/data  in  U.S.  Environmental 
Protection  Agency  (EPA)  reports  [EPA,  2004;  2012].  The  EPA’s  User's  Guide  for  Evaluating  Subsurface 
Vapor  Intrusion  into  Buildings  [EPA,  2004]  provides  insight  into  application  of  the  Johnson  and  Ettinger 
model,  including  relevant  calculations,  appropriate  ranges  of  parameters,  and  selected  default  values. 
Table  5  lists  relevant  ranges  and  default  values  from  EPA  [2004].  The  2012  Vapor  Intrusion  Database 
[EPA,  2012]  provides  a  data  set  from  which  inferences  can  be  drawn  regarding  “typical”  values  observed 
for  various  vapor  intrusion  related  parameters. 

Table  5.  Typical  ranges  of  parameter  values  for  residential  buildings  taken  from  EPA  [2004]  guidance 


Parameter 

Units 

Nominal  Range 
EPA  2004 

EPA  Recommended 
Default  Value 

building  footprint 

m2 

80  -  225 

100 

mixing  height 

m 

2.13-4.88 

2.44  (slab) 

3.66  (basement) 

air  exchange  rate 

day1 

2.4-36 

6 

h  (—  Ac  /  Ab) 

— 

0.0001  -0.0019 

0.000344  (slab) 
0.00022  (basement) 

Wc 

m 

0.0005-0.01 

0.001 

Os 

m3/d 

1.44-14.4 

7.2 

kv 

m2 

1.0  E-10-  1.0  E-1 6 

1.0  E-1 2 
(fine  sand) 

AP 

Pa 

0-20 

4 

Lc 

m 

0.1 

dz 

m 

0.15  (slab) 

2.0  (basement) 

EPA  [2004]  guidance  for  the  foundation/wall  thickness  simply  recommends  a  value  of  0. 1  m.  However, 
for  flexibility,  the  VIETUS  software  allows  values  up  to  0.61  m. 

The  building  ventilation  rate,  Qb ,  is  a  user  input  that  can  be  estimated  (externally  by  the  user)  from  other 
quantities  using  Equation  10.  Based  on  EPA  [2004]  information  (also  listed  in  Table  5),  a  typical 
residential  building  would  have  a  ventilation  rate  in  the  range  of  1500  to  2400  m3/d.  The  VIETUS 
software  suggests  a  range  of  1200  to  10,000  m3/d,  but  allows  the  user  to  specify  a  value  up  to  80,000  m3/d 
to  allow  flexibility  for  different  building  configurations. 

Qh  =  (footprint  area)  •  ( mixing  height)  ■  (air  exchange  rate)  ( 10) 
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The  A  inlet  term  of  Equation  9  represents  the  area  through  which  vapor  can  migrate  into  the  building.  For 
the  porous  foundation  case,  this  area  is  the  sub-grade  surface  area  of  the  foundations,  Ab ,  calculated  as 
shown  in  Equation  1 1 .  For  the  case  of  foundation  cracks,  the  area  of  the  cracks,  Ac ,  can  be  estimated 
either  from  the  crack  ratio,  q,  or  from  a  nominal  crack  width,  wc,  and  a  cumulative  crack  length,  xc,  as 
shown  in  Equations  12  and  13,  respectively  (and  where  rc  is  the  nominal  crack  radius).  A  common 
assumption  is  that  the  primary  cracks  are  at  the  joints  between  the  foundation  floor  and  the  walls.  Such  a 
“perimeter  crack”  would  have  a  total  length  (xc)  equal  to  the  building  perimeter,  as  calculated  from  the 
building  width  (Wb)  and  length  (Lb)  in  Equation  14.  The  calculations  in  Equations  1 1  to  14  are  used  in 
the  VIETUS  software  to  provide  defaults  based  on  the  secondary  input  parameter  values,  though  the  user 
can  override  the  default  calculations  by  entering  a  site-specific  number.  Permissible  ranges  are  not 
applied  directly  to  Ainiet.  Ab,  or  Ac,  though  constraints  on  q,  Wb,  Lb,  and/or  wc  will  constrain  the 
calculated  values  for  the  area.  In  addition,  a  limit  related  to  calculation  of  the  Qs  value  is  imposed  on  Ac , 
as  discussed  below. 


Ah  =  dx  ■  dy  +  2  ■  dx  ■  dz  +  2  ■  dy  ■  dz 

(11) 

Sr 

II 

4 

(12) 

4  =  wc-xc=2-rc-xc 

(13) 

xc=2-Wb+2-Lb 

(14) 

The  gas  flow  rate  into  the  building  from  the  surrounding  soil  can  either  be  directly  input  or  can  be 
estimated  using  the  approach  described  by  Johnson  and  Ettinger  [1991]  in  their  sample  calculations  for 
the  scenario  of  vapor  intrusion  through  cracks  in  the  foundation.  Johnson  and  Ettinger  assumed  that  the 
cracks  were  filled  with  dust/dirt  that  was  similar  in  nature  (e.g.,  porosity,  moisture  content)  to  the  soil 
underlying  the  foundation.  As  a  simplification,  the  approach  proposed  by  Nazaroff  [1988]  was  used  to 
estimate  Qs  from  a  total  crack  length  (e.g.,  representing  the  total  perimeter  for  the  floor/wall  seam),  the 
crack  width,  soil  properties,  and  building  under-pressurization.  This  approach  to  estimating  gas  flow  into 
the  building  is  based  on  an  idealized  model  of  cracks  at  the  floor/wall  seams,  where  the  crack  is  a  long 
cylinder  of  specified  radius.  Taking  the  crack  radius,  rc,  to  be  half  of  the  quotient  of  the  area  of  the 
cracks,  Ac ,  divided  by  the  cumulative  crack  length,  xc,  the  estimate  of  Qs  is  determined  from  Equation  15 
(a  variant  of  Johnson  and  Ettinger’ s  equation  24).  The  jjK  parameter  is  the  viscosity  of  air  at  the 
subsurface  temperature,  n  has  the  standard  meaning  of  the  constant,  and  the  other  parameters  for  this 
equation  are  identified  in  Table  2.  This  method  of  estimating  Qs  is  presented  as  an  optional  approach  that 
may  be  selected  if  the  assumptions  are  suitable  for  a  given  site,  or  which  can  be  neglected  in  favor  of 
direct  entry  of  Qs  (e.g.,  for  a  porous  foundation  scenario  or  based  on  measurements  from  tracer  tests). 

=  2-/r-AP-kv _ A, _ 

jug  In  (4-dz-xc)-\n(Ac)  ('15" 

As  with  the  Qb  calculation,  information  from  EPA  documents  is  used  as  a  guide  for  appropriate  ranges  of 
secondary  input  parameters.  The  air  viscosity  is  a  simple  lookup/interpolation  calculation  based  on 
published  data  [Lemmon,  2014].  The  permeability  of  soil  adjacent  to  the  foundation  is  very  site-specific, 
but  EPA  [2004]  provides  a  range  (see  Table  5)  and  recommends  a  default  value  of  1.0E-12  m2, 
representing  fine  sand.  Typical  ranges  of  building  under-pressure  are  also  provided  by  EPA  [2004] 
(Table  5).  The  EPA  guidance  is  also  used  to  constrain  Qs  values  themselves.  EPA  [2004]  suggests  that  a 
range  from  1.44  to  14.4  m3/d  is  reasonable,  while  a  review  of  data  [EPA,  2012]  indicates  that  Qs  is 
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typically  3%  or  less  of  the  building  ventilation  rate,  Qb .  The  VIETUS  software  suggests  an  upper  limit  of 
(xc/40)-15  to  allow  scaling  to  larger  building  sizes.  However,  Qs  values  as  large  as  6%  of  Qb  are 
permitted  (which  is  twice  the  observed  ratio  of  Qs/Qb). 


The  final  parameter  used  in  Equation  9  (and  not  listed  in  Table  2)  for  the  vapor  intrusion  estimate  is  the 
gas  diffusion  coefficient,  Dc ,  describing  vapor  diffusion  of  the  contaminant  through  the  porous  media  of 
the  cracks  or  building  foundation.  Using  the  same  Millington-Quirk  approach  as  taken  by  Johnson  and 
Ettinger  [1991],  the  value  of  Dc  can  be  calculated  from  Equation  16.  Here,  6c,g  is  the  gas  porosity  of  the 
cracks/foundation,  Qc, total  is  the  total  porosity  within  the  cracks/foundation,  and  Da  is  the  diffusion 
coefficient  in  air  (m2/day)  for  the  contaminant  of  interest.  The  EPA  [2004]  guidance  takes  the  crack 
moisture  content  to  be  “dry”  though  Johnson  and  Ettinger  [1991]  used  a  value  similar  to  the  underlying 
soil  in  their  example.  However,  the  Johnson  and  Ettinger  model  only  accounts  for  gas  diffusion  and 
convective  vapor  movement  into  a  building,  so  the  VIETUS  software  will  flag  the  6c,g  if  it  is  less  than 
50%  of  the  total  crack/foundation  porosity. 


d  ■ e 3 

£)  _  a  c,g 


f)2 

total 


(16) 


The  diffusion  coefficient  in  air,  which  is  used  to  determine  a  gas  diffusion  coefficient  for  the  porous 
cracks/foundation,  can  be  calculated  from  Equation  17  using  the  tabulated  contaminant-specific 
correlation  coefficients  listed  in  Table  6  [Yaws,  2003]. 

Da  =(A  +  B-T  +  C-T2)- (86400  /1002)  TisinK  (17) 


Table  6.  Tabulated  correlation  coefficients  for  calculating  the  gas  diffusion  coefficient  for  contaminants 
of  interest,  from  Yaws  [2003] 


Correl.  Coef.  for  Calculating  Gas  Diffusion  Coef. 

Contaminant 

Abbrev. 

Contaminant 

A 

B 

C 

CT 

Carbon  Tetrachloride 

-4.1620  E-2 

2.5749  E-4 

5.3598  E-7 

CF 

Chloroform 

-4.9960  E-2 

2.8829  E-4 

5.9407  E-7 

DCM 

Dichloromethane 

-5.6860  E-2 

3.3534  E-4 

6.8151  E-7 

CM 

Chloromethane 

-6.9250  E-2 

4.2133  E-4 

8.2225  E-7 

PCE 

Tetrachloroethene 

-3.6840  E-2 

2.3847  E-4 

5.1123  E-7 

TCE 

Trichloroethene 

-4.1270  E-2 

2.6580  E-4 

5.5704  E-7 

1,1-DCE 

1,1-Dichloroethene 

-4.5030  E-2 

2.9721  E-4 

6.0085  E-7 

cDCE 

cis-1 ,2-Dichloroethene 

-4.6220  E-2 

3.0091  E-4 

6.1981  E-7 

tDCE 

trans-1 ,2-Dichloroethene 

-4.6380  E-2 

3.0373  E-4 

6.2044  E-7 

VC 

Vinyl  Chloride 

-5.4240  E-2 

3.6658  E-4 

7.2028  E-7 

1 ,1 ,1 ,2-TeCA 

1,1,1 ,2-T  etrachloroethane 

-3.6290  E-2 

2.2950  E-4 

4.9507  E-7 

1 ,1 ,2,2-TeCA 

1 ,1 ,2,2-Tetrachloroethane 

-4.0420  E-2 

2.2899  E-4 

4.9882  E-7 

1,1,1-TCA 

1,1,1-Trichloroethane 

-4.5040  E-2 

2.5770  E-4 

5.3553  E-7 

1,1,2-TCA 

1 ,1 ,2-Trichloroethane 

-4.4400  E-2 

2.5339  E-4 

5.4055  E-7 

1,1-DCA 

1,1-Dichloroethane 

-4.7450  E-2 

2.8979  E-4 

5.9566  E-7 

1,2-DCA 

1,2-Dichloroethane 

-4.8920  E-2 

2.8847  E-4 

6.0315  E-7 

CE 

Chloroethane 

-5.5930  E-2 

3.4520  E-4 

6.8934  E-7 

Equations  9  and  15  both  include  terms  that  can  cause  the  calculations  to  fail.  Thus,  the  set  of  parameters 
that  comprise  these  terms  are  constrained  to  avoid  calculation  failure.  For  Equation  15,  the  ratio 
4-dz-Xc  /  Ac  must  be  greater  than  1.  The  user  can  avoid  this  calculational  constraint  by  directly  entering  a 
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Qs  value  (setting  the  secondary  parameters  xc  and  Ac  to  empty  values).  For  Equation  9,  the 
(Qs-Lc)/(Dc-  Ainiet)  term  must  be  of  a  magnitude  to  avoid  failure  of  the  exponential  function.  If  the  site- 
specific  configuration,  as  entered  into  the  VIETUS  software,  results  in  a  failure  due  to  this  term  (i.e.,  for 
convection  dominated  transport  through  the  foundation  cracks),  then  the  user  is  referred  to  Johnson  and 
Ettinger  [1991],  who  discuss  an  alternate  formulation  for  their  attenuation  factor  (in  their  equation  22). 

2.5  Limitations  and  Assumptions 

The  procedure  described  in  this  report  provides  an  estimate  of  vapor  intrusion  based  on  user-defined 
inputs  and  a  set  of  constraints/assumptions.  The  usefulness  of  such  estimates  must  be  evaluated  for  a 
specific  site,  considering  whether  the  constraints/assumptions  are  appropriate.  One  element  to  consider  is 
the  uncertainty  of  input  parameter  values.  While  input  parameters  are  often  not  known  precisely,  a 
sensitivity  analysis  can  readily  be  conducted  with  the  VIETUS  software  to  assess  the  potential  impact  of 
reasonable  variation  in  specific  input  parameters  on  the  estimated  vapor  intrusion  results.  Another  aspect 
to  consider  is  the  relevance  of  factors  outside  the  scope  of  the  calculational  approach  used  in  the  VIETUS 
software,  which  constrains  the  scenarios  for  which  the  estimates  are  suitable.  Other  evaluation 
approaches  (perhaps  in  conjunction  with  estimates  from  the  VIETUS  software)  may  be  needed  to 
consider  factors  such  as  the  degree  of  source  depletion  over  time,  adsorption,  biological  transformation, 
other  physical  attenuation  mechanisms,  multiple  sources,  sources  in  the  groundwater,  and  variation  in 
recharge/infiltration  through  the  vadose  zone  source.  In  addition  to  constraints  on  factors  include  in  the 
calculational  approach,  the  appropriateness  of  simplifying  assumptions  used  in  the  approach  should  be 
considered  with  respect  to  the  site-specific  conditions.  For  instance,  the  generalized  conceptual  model 
used  in  the  approach  is  appropriate  for  sites  where  vapor-phase  transport  dominates  contaminant 
movement. 

It  is  important  to  understand  the  scope  of  the  calculation  approach/tool  presented  here.  Figure  1  depicts 
the  vapor-phase  diffusion  from  the  current  vadose  zone  source  and  the  resulting  steady-state  distribution 
of  contamination  in  the  vadose  zone  and  the  associated  contamination  in  groundwater.  VIETUS  is  not 
configured  to  address  vapors  that  arise  either  solely  from  groundwater  plumes  (i.e.,  where  contamination 
was  released  directly  to  groundwater  or  where  a  vadose  zone  source  is  not  currently  present)  or  from 
groundwater  contamination  resulting  from  a  historically  stronger  vadose  zone  source.  This  tool  is 
intended  to  supply  information  to  support  a  decision  about  SVE  disposition  relative  to  diminishing  returns 
in  remediation  of  a  vadose  zone  contaminant  source  (i.e.,  “If  we  stop  SVE  now,  would  we  expect  a  vapor 
intrusion  issue  from  the  remaining  source?”). 

Specific  Limitations  and  Assumptions: 

•  Vapor-phase  transport  dominates  contaminant  movement  in  the  vadose  zone  source. 

•  The  specific  site  can  be  suitably  represented  by  the  generalized  conceptual  site  model. 

•  The  contaminant  source  can  be  represented  as  single  source  area  in  the  vadose  zone  within  a 
square  footprint  of  specified  size  and  thickness.  Note  that  it  may  be  appropriate  for  a  cluster  of 
separate  sources  to  be  represented  as  a  single  source  with  respect  to  vapor  transport  (see  Truex  et 
al.  [2013]  and  the  discussion  of  the  generalized  conceptual  model  in  Section  2.0). 

•  The  contaminant  source  strength  can  be  represented  as  a  constant  vapor-phase  contaminant 
concentration  or  a  mass  discharge  of  contaminant  that  results  in  a  steady-state  distribution  of 
contaminant  in  the  subsurface. 
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•  Simulations  were  run  to  obtain  steady-state  concentration  distributions  in  the  gas  and  aqueous 
phases  throughout  the  computational  domain. 

•  There  are  no  additional  sources  in  the  vadose  zone  or  groundwater. 

•  The  subsurface  can  be  approximated  as  a  homogeneous  material  (i.e.,  is  not  layered  or 
heterogeneous)  with  uniform  properties  (temperature,  vadose  zone  moisture  content,  bulk  density, 
porosity).  Note  that  some  types  of  heterogeneity  may  have  minimal  impact  on  vapor  transport 
and  can  be  approximated  as  homogeneous  (see  Truex  et  al.  [2013]  and  the  discussion  of  the 
generalized  conceptual  model  in  Section  2.0). 

•  The  transverse  transport  of  contaminant  (perpendicular  to  the  numerical  model  x-axis)  is 
symmetrical,  so  only  half  of  a  model  domain  in  the  y-axis  direction  need  be  modeled.  Note  that 
the  direction  of  groundwater  flow  is  parallel  to  the  numerical  model  x-axis. 

•  Extrapolation  outside  the  bounds  of  the  model  domain  or  bounds  of  parameter  key  values  is  not 
allowed.  See  Table  1  and  Table  2  for  lists  of  the  permissible  ranges  for  parameters. 

•  Approximations  inherent  in  using  linear  interpolation  between  key  values  are  small  relative  to  the 
accuracy  needed  for  the  estimate. 

•  Differences  in  contaminant-specific  gas  phase  diffusion  coefficients  are  small  enough  to  neglect 
when  scaling  the  soil-gas  concentration  results  from  the  pre -modeled  scenarios. 

•  Soil  gas  concentration  data  from  the  simulations  consists  of  values  greater  than  1  ng/L  (values 
below  that  cutoff  are  ignored). 

•  Foundations  do  not  extend  more  than  4  m  below  ground  surface  and  Vi  m  depth  increments 
provide  sufficient  resolution  for  specifying  a  site  scenario. 

•  Recharge  is  assumed  to  be  low,  such  that  use  of  simulation  results  based  on  a  recharge  of  0.4 
cm/yr  is  appropriate  to  represent  conditions  where  vapor-phase  transport  in  the  vadose  zone 
dominates  (see  also  the  discussion  below). 

•  Vapor  intrusion  is  based  on  an  adaptation  of  the  Johnson  and  Ettinger  [1991]  model.  Numerical 
simulation  results  are  used  for  the  subsurface  transport  component  while  the  Johnson  and  Ettinger 
model  is  applied  to  determine  the  entry  rate  of  contaminant  into  the  building  (as  driven  by 
building  ventilation)  and  the  associated  building  concentration  for  the  contaminant.  This 
decoupling  of  diffusive  flux  in  the  subsurface  and  the  rate  of  vapor  intrusion  will  result  in  vapor 
intrusion  estimates  that  differ  from  the  full  Johnson  and  Ettinger  model. 

•  The  distribution  of  soil-gas  concentration  in  the  subsurface  is  a  result  of  diffusive  transport  and 
does  not  account  for  advection  (from  barometric  pumping,  density  effects,  or  induced  by  the 
building  under-pressurization.  However,  the  EPA  [2004]  note  that  the  zone  of  influence  for 
advective  contaminant  flow  due  to  vapor  intrusion  is  expected  to  be  limited  to  a  region 
immediately  adjacent  to  the  building  foundation. 

•  A  single  soil-gas  contaminant  concentration  is  used  to  determine  the  concentration  in  the 
building. 

•  Vapor  intrusion  occurs  through  either  cracks  (which  could  be  floor/wall  seams)  or  a  porous 
foundation. 
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•  Cracks  are  assumed  to  be  filled  with  dust/dirt  with  sufficient  air-filled  porosity  to  allow  vapor 
intrusion  into  the  building. 

•  The  estimate  of  gas  flow  into  the  building  based  on  secondary  vapor  intrusion  parameters 
assumes  an  idealized  model  of  cracks  at  the  floor/wall  seams,  where  the  crack  is  a  long  cylinder 
of  specified  radius  [Johnson  and  Ettinger,  1991].  However,  the  user  can  externally  determine  this 
gas  flow  into  the  building  and  can  input  a  specific  value. 

•  Gas  flow  into  a  building  for  a  porous  foundation  occurs  through  the  entire  subsurface  area  of  the 
foundation  [Johnson  and  Ettinger,  1991]. 

•  The  temperature  dependence  of  the  Henry’ s  Law  constant  is  determined  as  suggested  by  Brennan 
et  al.  [1998]  using  vapor  pressure  and  solubility  values. 

•  The  Antoine  correlation  [Yaws  et  al.,  2009]  is  suitable  for  determining  vapor  pressure  at  a 
specified  temperature. 

•  The  correlations  of  Mackay  et  al.  [2006]  and  Yaws  [2003]  are  suitable  for  determining  solubility 
at  a  specified  temperature. 

•  The  correlations  of  Yaws  [2003]  are  suitable  for  determining  the  gas  diffusion  at  a  specified 
temperature. 

The  impact  of  recharge  on  the  soil-gas  contaminant  concentration  in  the  vadose  zone  is  likely  to  be 
complicated.  Recharge  affects  the  mass  flux  from  the  vadose  zone  source  into  the  groundwater.  Once  in 
the  groundwater,  contaminant  mass  can  be  transported  downgradient  and  then  volatilize  back  into  the 
vadose  zone  at  a  distance  from  the  source.  This  effect  could  occur  due  to  higher  recharge,  a  strong  source 
close  to  the  groundwater,  or  a  source  within  the  groundwater.  Regardless,  the  relationships  are 
complicated  and  it  is  possible  that  a  significant  amount  of  additional  simulations  would  be  required  to 
understand  the  behavior.  Thus,  the  VIETUS  software  currently  focuses  only  on  estimating  the  impact  of 
the  vadose  zone  source  on  vapor  intrusion  for  a  low  value  of  recharge  where  vapor-phase  transport  in  the 
vadose  zone  dominates.  Site-specific  simulations  should  be  investigated  for  cases  where  recharge  is  more 
significant  (or  where  a  significant  source  within  the  groundwater  exists). 
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3.0  Example  Calculation 


To  illustrate  the  procedure  for  estimating  the  soil  gas  and  building  contaminant  concentrations,  this 
example  walks  through  the  calculations  for  the  scenarios  shown  in  Table  7.  These  scenarios  represent  the 
compiled  set  of  input  data  for  three  variants  (Cases  A,  B,  and  C),  which  differ  in  the  way  that  the  source 
strength  is  specified,  the  way  in  which  the  crack  area  is  calculated,  the  method  for  input  of  Qs ,  and  the 
lateral  location  of  the  building.  Case  A  uses  a  crack  width  to  estimate  the  area  of  cracks  and  has  a  directly 
specified  value  of  Qs .  Case  B  uses  a  crack  ratio  to  estimate  the  area  of  cracks  and  Qs  is  from  the  default 
calculation.  Case  C  is  similar  to  Case  B,  but  the  source  strength  is  specified  as  a  mass  discharge,  the 
lateral  location  of  the  building  is  closer  to  the  source,  and  the  under-pressure  is  halved. 


Table  7.  User  input  for  the  scenario  variants  applied  in  the  example  calculations  for  estimation  of  the 
soil  gas  concentrations  and  vapor  intrusion 


User  Input  -  Source/Transport  Parameters 

Scenario  Name: 

— 

Case  A 

Case  B 

Case  C 

Contaminant: 

— 

TCE 

TCE 

TCE 

T 

Temperature: 

ra 

20 

20 

20 

to 

Avg.  Moisture  Content: 

[wt  %] 

2.825 

2.825 

2.825 

VZT 

Vadose  Zone  Thickness: 

[m] 

20 

20 

20 

LI 

Depth  to  Top  of  Source: 

[m] 

5.85 

5.85 

5.85 

Z 

Source  Thickness: 

[ml 

3.5 

3.5 

3.5 

w(=l) 

Source  Width  (=  Length): 

[m] 

15.8 

15.8 

15.8 

q 

GW  Darcy  Velocity: 

fm/davl 

0.165 

0.165 

0.165 

Source  Strength  Input  Type: 

— 

Gas  Concentration 

Gas  Concentration 

Mass  Discharge 

Cqs 

Source  Gas  Concentration: 

[ppmv] 

159 

159 

Msrc 

Source  Mass  Discharge: 

[g/day]  ■■■■ 

Total  porosity  =0.3.  dry  bulk  density  =  1.855  g/mL,  and  recharge  =0.4  cmfyr 

User  Input  -  Vapor  Intrusion  Parameters 

dx 

Longitudinal  Distance  to  Building: 

[m] 

50 

50 

20 

dy 

Transverse  Distance  to  Building: 

[m] 

10 

10 

20 

dz 

Depth  of  Basement/Foundation: 

M 

4 

4 

4 

Lc 

Foundation/Wall  Thickness: 

[m] 

0.1 

0.1 

0.1 

Q0 

Building  Ventilation  Rate: 

[rrf/day] 

2000 

2000 

2000 

User  Input  —  Vapor  Intrusion  Parameters  (directly specified  or  calculated  from  2°  parameters) 

Qs 

Gas  Flow  Into  Building  from  Soil: 

[nf/day] 

8.0 

1.862 

0.931 

Anlet 

Inlet  Area  (Cracks  or  "Pores"): 

[m2] 

0.6 

0.44 

0.44 

9c.  total 

Total  Porosity  within  Crack/Fndn: 

[~] 

0.3 

0.3 

0.3 

®c.g 

Gas-Filled  Porosity  for  Cracks/Fndn: 

[-] 

0.225 

0.225 

0.225 

User  Input- Secondary  Vapor  Intrusion  Parameters  (forderiving  Q,  and  An«) 

Vapor  Intrusion  Path  Type: 

— 

Cracks 

Cracks 

Cracks 

w„ 

Building  Width: 

[m] 

10 

10 

10 

Lb 

Building  Length: 

[m] 

20 

20 

20 

AP 

Building  Underpressure: 

[Pa] 

8 

4 

kv 

Near-Building  Soil  Gas  Permeability: 

[m2] 

1.00E-12 

1.00E-12 

At  from  Crack  Ratio/Width/Direct: 

— 

Width 

Ratio 

Ratio 

n 

Crack  Ratio  (Ac  /  Ab): 

H 

0.001 

0.001 

wc 

Crack  Width: 

[m] 

0.01 

*c 

Cumulative  Lateral  Crack  Length: 

[m] 

60 

60 

60 

Ac 

Total  Area  of  Foundation  Cracks: 

[m2] 

0.6 

0.44 

0.44 

The  calculation  procedure  begins  by  performing  the  interpolation,  which  is  done  in  two  phases:  a  spatial 
interpolation  phase  and  a  phase  to  interpolate  between  key- value  scenarios.  A  key- value  scenario  is 
defined  as  a  specific  permutation  of  key  values  for  the  six  source/transport  parameters  (RSP,  q,  SA,  VZT, 
STR7  and  S, ,  as  discussed  in  Table  1)  that  exhibit  a  nonlinear  concentration  response. 
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The  first  phase  interpolates  soil  gas  data  between  numerical  model  grid  cell  centers  and  model  layers  to 
obtain  a  soil  gas  concentration  for  the  arbitrary  lateral  location  ( dx  and  dy)  and  depth  of 
basement/foundation  ( dz )  that  are  specified  in  the  inputs.  Figure  6  depicts  both  a  set  of  concentrations 
from  the  simulation  results  (using  the  26th  key-value  scenario  of  parameters  for  Cases  A  and  B  as  an 
example)  and  the  process  of  interpolating  in  three  dimensions  to  obtain  the  concentration  at  the  point  of 
interest  (marked  with  the  star).  Concentrations  at  grid  cell  centers  (dots)  from  the  layer  above  and  the 
layer  below  the  depth  of  interest  are  interpolated  to  obtain  four  “corner”  concentrations  (squares)  at  the 
depth  of  interest.  Then  a  bilateral  interpolation  is  performed  to  determine  the  soil  gas  concentration,  Cs,m , 
at  the  lateral  location  of  interest  (star).  This  spatial  interpolation  process  is  repeated  to  supply  the  64  CSm 
values  used  in  the  second  phase  of  interpolation. 


867.8  273.0 


Figure  6.  Example  of  the  spatial  interpolation  from  numerical  grid  cell  centers  (dots)  to  the  location  of 
interest  (denoted  with  a  star)  within  the  bounds  of  those  cell  centers.  Interpolation  is  done 
vertically  first  to  obtain  the  four  corners  of  the  layer  of  interest  (squares),  then  done  laterally 
to  find  the  soil  gas  concentration  at  the  point  of  interest. 

The  second  phase  of  interpolation  involves  interpolating  between  key-value  scenarios  for  the  RSP,  q,  SA, 
VZT,  STR,  and  Sr  parameters  to  match  the  site-specific  conditions.  The  user  directly  specifies  the  site- 
specific  values  of  VZT  and  q,  while  values  of  RSP,  SA,  STR,  and  .S’,  (i.e.,  the  converted  value  of  co,  as 
discussed  in  the  footnotes  to  Table  1)  are  calculated  from  other  input  parameters.  Table  8  lists  the  site- 
specific  values  of  these  six  parameters  for  the  three  example  cases.  Each  of  the  64  key-value  scenarios 
between  which  the  interpolation  occurs,  represents  a  simulation  with  parameters  at  key  values  above  or 
below  (or  perhaps  equal  to)  the  user-specified  parameter  values.  Table  9  lists  the  CSim  values  for  Cases  A 
and  B  and  the  subsequent  interpolation  results  for  the  six  parameters  to  arrive  at  the  final  soil  gas 
concentration,  Cgu .  Case  C  has  different  values  because  of  the  different  spatial  location  specified. 

Table  8.  Parameter  values  for  the  example  cases  that  are  used  in  the  second  phase  of  interpolation 


Scenario  Name: 

— 

Case  A 

Case  B 

Case  C 

Sr 

Residual  Saturation: 

[--] 

0.175 

0.175 

0.175 

STR 

Source  Thickness  Ratio: 

[--] 

0.175 

0.175 

0.175 

VZT 

Vadose  Zone  Thickness: 

[m] 

20 

20 

20 

SA 

Areal  Footprint  of  Source: 

[m2] 

250 

250 

250 

q 

GW  Darcy  Velocity: 

fm/dayl 

0.165 

0.165 

0.165 

RSP 

Relative  Source  Position: 

[--] 

0.55 

0.55 

0.55 

24 


Table  9.  Sequential  interpolation  to  determine  Cgu  (ppbv)  for  Cases  A  and  B  of  the  example;  starting 
with  the  64  Csim  values,  results  of  linear  interpolations  between  adjacent  pairs  of  values  are 
shown  in  the  columns,  from  left  to  right,  with  the  rightmost  column  showing  the  value  of  Cgu 


Key- Value  Scenario 

Csim  ValUG 

After  RSP 

After  q 

After  SA 

After  VZT 

After  STR 

After  Sr 

(RSP ,  q. 

SA, 

VZT, 

STR, 

Sr) 

(ppbv) 

Interpolation 

Interpolation 

Interpolation 

Interpolation 

Interpolation 

Interpolation 

0.1,  0.03, 

100, 

10, 

0.1  , 

0.05 

28.0 

44.7 

1.0,  0.03, 

100, 

10, 

0.1  , 

0.05 

61.3 

64.2 

0.1,  0.3  , 

100, 

10, 

0.1  , 

0.05 

52.7 

83.7 

1.0,  0.3  , 

100, 

10, 

0.1  , 

0.05 

114.8 

128.8 

0.1,  0.03, 

400, 

10, 

0.1  , 

0.05 

97.2 

136.6 

1.0,  0.03, 

400, 

10, 

0.1  , 

0.05 

176.1 

193.8 

0.1,  0.3  , 

400, 

10, 

0.1  , 

0.05 

179.3 

251.0 

1.0,  0.3  , 

400, 

10, 

0.1  , 

0.05 

322.9 

485.6 

0.1,  0.03, 

100, 

30, 

0.1  , 

0.05 

339.5 

582.7 

1.0,  0.03, 

100, 

30, 

0.1  , 

0.05 

826.6 

580.4 

0.1,  0.3  , 

100, 

30, 

0.1  , 

0.05 

334.0 

578.1 

1.0,  0.3  , 

100, 

30, 

0.1  , 

0.05 

823.0 

842.4 

0.1,  0.03, 

400, 

30, 

0.1  , 

0.05 

775.9 

1109.1 

1.0,  0.03, 

400, 

30, 

0.1  , 

0.05 

1443.3 

1105.7 

0.1,  0.3  , 

400, 

30, 

0.1  , 

0.05 

763.3 

1102.4 

1.0,  0.3  , 

400, 

30, 

0.1  , 

0.05 

1442.5 

531.3 

0.1,  0.03, 

100, 

10, 

0.25, 

0.05 

22.2 

36.5 

1.0,  0.03, 

100, 

10, 

0.25, 

0.05 

50.8 

52.9 

0.1,  0.3  , 

100, 

10, 

0.25, 

0.05 

42.4 

69.4 

1.0,  0.3  , 

100, 

10, 

0.25, 

0.05 

96.5 

104.4 

0.1,  0.03, 

400, 

10, 

0.25, 

0.05 

70.7 

108.8 

1.0,  0.03, 

400, 

10, 

0.25, 

0.05 

147.0 

156.1 

0.1,  0.3  , 

400, 

10, 

0.25, 

0.05 

132.6 

203.3 

1.0,  0.3  , 

400, 

10, 

0.25, 

0.05 

274.3 

576.9 

0.1,  0.03, 

100, 

30, 

0.25, 

0.05 

443.8 

769.2 

1.0,  0.03, 

100, 

30, 

0.25, 

0.05 

1095.6 

766.8 

0.1,  0.3  , 

100, 

30, 

0.25, 

0.05 

436.4 

764.4 

1.0,  0.3  , 

100, 

30, 

0.25, 

0.05 

1093.4 

1049.5 

0.1,  0.03, 

400, 

30, 

0.25, 

0.05 

854.9 

1335.4 

1.0,  0.03, 

400, 

30, 

0.25, 

0.05 

1817.4 

1333.5 

0.1,  0.3  , 

400, 

30, 

0.25, 

0.05 

842.5 

1331.7 

1.0,  0.3  , 

400, 

30, 

0.25, 

0.05 

1822.3 

617.5 

0.1,  0.03, 

100, 

10, 

0.1  , 

0.3 

76.6 

120.9 

1.0,  0.03, 

100, 

10, 

0.1  , 

0.3 

165.4 

188.0 

0.1,  0.3  , 

100, 

10, 

0.1  , 

0.3 

160.3 

255.0 

1.0,  0.3  , 

100, 

10, 

0.1  , 

0.3 

350.0 

375.2 

0.1,  0.03, 

400, 

10, 

0.1  , 

0.3 

261.1 

363.9 

1.0,  0.03, 

400, 

10, 

0.1  , 

0.3 

467.1 

563.2 

0.1,  0.3  , 

400, 

10, 

0.1  , 

0.3 

542.1 

762.5 

1.0,  0.3  , 

400, 

10, 

0.1  , 

0.3 

983.7 

660.3 

0.1,  0.03, 

100, 

30, 

0.1  , 

0.3 

385.3 

670.0 

1.0,  0.03, 

100, 

30, 

0.1  , 

0.3 

955.7 

652.6 

0.1,  0.3  , 

100, 

30, 

0.1  , 

0.3 

358.0 

635.2 

1.0,  0.3  , 

100, 

30, 

0.1  , 

0.3 

913.4 

945.4 

0.1,  0.03, 

400, 

30, 

0.1  , 

0.3 

877.4 

1268.0 

1.0,  0.03, 

400, 

30, 

0.1  , 

0.3 

1659.8 

1239.6 

0.1,  0.3  , 

400, 

30, 

0.1  , 

0.3 

817.0 

1211.1 

1.0,  0.3  , 

400, 

30, 

0.1  , 

0.3 

1606.5 

703.8 

0.1,  0.03, 

100, 

10, 

0.25, 

0.3 

60.7 

99.2 

1.0,  0.03, 

100, 

10, 

0.25, 

0.3 

137.8 

155.9 

0.1,  0.3  , 

100, 

10, 

0.25, 

0.3 

128.8 

212.6 

1.0,  0.3  , 

100, 

10, 

0.25, 

0.3 

296.7 

306.7 

0.1,  0.03, 

400, 

10, 

0.25, 

0.3 

191.1 

292.6 

1.0,  0.03, 

400, 

10, 

0.25, 

0.3 

394.5 

458.3 

0.1,  0.3  , 

400, 

10, 

0.25, 

0.3 

403.4 

623.9 

1.0,  0.3  , 

400, 

10, 

0.25, 

0.3 

845.0 

747.3 

0.1,  0.03, 

100, 

30, 

0.25, 

0.3 

503.3 

888.2 

1.0,  0.03, 

100, 

30, 

0.25, 

0.3 

1274.3 

869.2 

0.1,  0.3  , 

100, 

30, 

0.25, 

0.3 

469.7 

850.2 

1.0,  0.3  , 

100, 

30, 

0.25, 

0.3 

1231.9 

1187.9 

0.1,  0.03, 

400, 

30, 

0.25, 

0.3 

966.2 

1531.7 

1.0,  0.03, 

400, 

30, 

0.25, 

0.3 

2099.0 

1508.0 

0.1,  0.3  , 

400, 

30, 

0.25, 

0.3 

907.9 

1484.3 

1.0,  0.3  , 

400, 

30, 

0.25, 

0.3 

2062.5 
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The  first  column  of  Table  9  lists  the  C«m  values  from  the  spatial  interpolation  for  Cases  A  and  B.  Each 
subsequent  column  is  determined  by  applying  the  linear  interpolation  of  Equation  1  to  adjacent  values. 
The  key  RSP  values  (Table  1)  are  0.1  and  1.0,  because  the  site-specific  RSP  value  is  0.55  (Table  8).  The 
first  pair  of  Cs/m  values  in  Table  9  corresponds  to  those  two  key  values  of  RSP  (for  the  same  key  values  of 
q,  SA,  VZT ',  STR,  and  Sr),  so  interpolation  for  RSP  can  be  applied  as: 


0.55-0.1 

1.0-0.1 


■(61. 3 -28.0) +  28.0 


44.7  ppbv 


This  same  interpolation  process  is  performed  for  the  remaining  pairs  of  CSm  values  to  get  the  intermediate 
interpolation  results  in  the  second  column  of  Table  9.  Pairs  of  values  from  the  second  column  are  then 
interpolated  to  obtain  concentration  values  for  the  user-specified  q  value  of  0.165  m/day  (Table  8).  The 
key  q  values  (Table  1)  above  and  below  the  site-specific  q  are  0.03  and  0.3  m/day.  The  interpolation  of 
the  first  pair  of  concentrations  in  the  second  column  thus  gives: 


0.165-0.03 
0.3  -  0.03 


•  ( 44.7  -  83.7)  +  83.7  =  64.2  ppbv 


This  process  continues  in  the  same  fashion  for  pairs  of  values  with  interpolation  between  key  values  to 
the  site-specific  conditions  for  the  q,  SA ,  VZT.  STR ,  and  S,  parameters.  The  results  of  these  sequential 
linear  interpolations  are  shown  in  Table  9,  with  the  final  interpolation  providing  the  Cgu  value  of  617.5 
ppbv  (for  Cases  A  and  B). 

For  Case  C  of  the  example,  interpolation  is  required  to  determine  the  source  mass  discharge,  MSim, 
corresponding  to  the  site-specific  parameter  values.  This  value  is  used  in  the  scaling  process. 
Interpolation  to  determine  MSim  does  not  involve  the  first  phase  of  spatial  interpolation  (because  it  pertains 
to  the  defined  source),  but  does  involve  the  second  phase  to  interpolate  to  the  site-specific  conditions  for 
the  RSP,  q,  SA,  VZT,  STR,  and  Sr  parameters  in  a  process  parallel  to  that  shown  in  Table  9.  The  initial 
values  of  mass  discharge,  from  which  the  64  values  are  drawn,  are  tabulated  in  the  VIETUS  software  (see 
the  SVEET  documentation  [Truex  et  al.,  2013]  for  a  list  of  all  mass  discharge  values).  For  Case  C,  the 
resultant  Msm  is  58.3  g/day  and  the  Cgu  value  is  1425.9  ppbv. 

Once  Cgi,  (and  M,,m ,  if  needed)  are  determined,  the  next  step  is  to  apply  scaling  factors  for  the  two 
parameters  (Henry’s  Law  Constant  and  source  strength)  where  concentration  exhibits  a  linear  response. 

The  vapor  pressure  and  solubility  for  trichloroethene  at  20  °C  are  calculated  using  contaminant-specific 
correlations  in  Equations  2  and  3  (with  conversion  of  solubility  to  mole  fraction  based  on  Equation  5): 

Pvap  =  10[6-87981  -  1 157.83  /  (20.0  +  202.58)]  =  47.6376  mm  Hg 

xp  =  [1.4049  -  0.0082223- (20.0  +  273.15)  +  0.000013218-(20.0  +  273.15)2]  /  100  =  0.13045 
x  =  (0.13045  /  100M18.01528/  131.387)=  1.7886xl04 

These  results  are  used  with  Equation  6  to  obtain  the  dimensionless  value  for  the  Henry’s  Law  constant: 

47.6376  18.01528  1  1 

- =  0.263 

0.00017886  0.9982-0.08205746-293.15  1000  760 
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For  Cases  A  and  B  of  the  example,  Equation  7  is  applied  to  obtain  the  final  estimated  soil  gas 
contaminant  concentration,  Cg ,  at  the  specified  location  beneath  the  building  foundation.  For  Case  C, 
Equation  8  is  applied  to  obtain  Cg  because  the  source  strength  is  specified  as  a  mass  discharge.  These 
final  calculations  give  the  following  results: 

Cases  A&B :  617.5  ■  ■  ^9-®  _  2090  ppbv 

0.263  159.0 

Case  C  :  1425.9  ■  °^9-  ■  —  =  4798  ppbv 

0.263  58.3 


Given  the  soil  gas  concentration  beneath  the  building  foundations,  the  final  step  is  to  generate  the 
estimate  of  the  contaminant  concentration  in  the  building  due  to  vapor  intrusion.  The  diffusion 
coefficient  in  air  for  TCE  at  20  °C  is  calculated  with  the  correlation  formula  of  Equation  17  and 
associated  constants  in  Table  6: 

Da  =  -4. 127E-2  +  (2.658E-4)(293.15)  +  (5.5704E-7)(293.15)2  =  0.0845  cm2/s  =  0.73025  m2/day 


The  diffusion  coefficient  through  cracks  in  the  foundation  is  then  calculated  with  Equation  16: 


D. 


10. 

0.73025  -(0.225) 3 

oJ2 


0.0562  m2  /  day 


The  area  through  which  vapor  intrusion  occurs  64  „,/,.,)  is  the  area  of  the  cracks  64,.)  for  the  three  cases 
considered  here.  For  Case  A,  a  crack  width  is  specified  and  the  crack  area  is  calculated  as  the  crack  width 
times  the  cumulative  lateral  crack  length  (i.e.,  the  building  perimeter).  For  Cases  B  and  C,  the  crack  area 
is  estimated  based  on  a  crack  ratio  and  the  total  surface  area  of  the  below-ground  walls  and  floor. 

Case  A:  Ainie,  -  Ac  -  wc-xc  =  wc-(2-  Wb  +  2-Lb )  =  (0.01)  (60)  =  0.6  m2 

Cases  B  &  C:  Am/e/  -  Ac  =  rj-Ab  =  tj-( 2-  Wb-dz  +  2 -Lb-dz  +  Wb-Lb )  =  (0.001)  (440)  =  0.44  m2 


Case  A  uses  a  directly  specified  value  for  the  gas  flow  into  the  building.  Cases  B  and  C  use  Equation  15 
to  estimate  Qs .  For  Case  B  the  calculation  is: 


2-71-8  ■  10'12 
2.107x10  10 


_ 60 _ 

ln(4  •  4  ■  60)  -  \n(0.44) 


=  1.862  m  3 1  day 


The  same  calculation  with  Case  C  parameter  values  gives  Qs  =  0.931  (because  the  building  under¬ 
pressure  is  half  as  much). 


Having  defined  all  key  values,  Equation  9  can  be  applied  to  determine  the  estimate  for  the  contaminant 
concentration  in  the  building.  First  calculating  the  term  inside  the  exponential  function: 

QpL 

— — -  =  [(0.1)  (8)]/ [(0.6)  (0.0562)]  =  23.72  forCaseA 

Dc-A 

=  [(0.1)  (1.862)]/ [(0.44)  (0.0562)]  =  7.53  forCaseB 

=  [(0.1)  (0.931)]/ [(0.44)  (0.0562)]  =  3.76  for  Case  C 
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Then,  completing  the  calculation  for  each  separate  example  case,  the  contaminant  concentration  in  the 
building  is  estimated  to  be: 


8  ■  2090  •  exp(23.72) 


Ch  2000  ■  exp(2J.72)  -  2000  +  8 


=  8.4  ppbv 


1.862  ■  2090  •  exp(7.5i) 


C  = _ 

b  2000  ■  exp(7.5T)  -  2000  + 1 .862 


0.931-  4798  -exp(i. 76) 


Q  - 

h  2000  ■  exp(3.76)~  2000  +  0.931 


=  2.0  ppbv 


=  2.3  ppbv 


for  Case  A 


for  Case  B 


for  Case  C 
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4.0  VIETUS  User  Guide 


For  user  convenience,  the  calculational  procedure  for  estimating  a  soil-gas  contaminant  concentration  in 
the  vadose  zone  and  the  associated  vapor  intrusion  concentration  (described  in  Section  2.0  and  illustrated 
in  Section  3.0)  has  been  implemented  in  a  spreadsheet  software  tool.  The  VIETUS  software  allows  the 
user  to  easily  enter  data  and  calculate  the  estimated  soil  gas  and  building  concentrations  for  one  or  more 
scenarios  conforming  to  the  generalized  conceptual  model  described  in  Section  2.0.  The  system 
requirements,  installation,  user  interface,  and  application  of  the  VIETUS  software  are  described  in  the 
sections  below. 

4.1  System  Requirements 

The  following  hardware  and  software  are  required  to  use  the  VIETUS  software: 

•  Personal  computer  based  on  Intel®  IA-32  or  Intel®  64  processor  architectures, 

•  Microsoft®  Windows®  XP  or  Microsoft®  Windows®  7  operating  system, 

•  Microsoft®  Excel®  2007  or  Excel®  2010 

The  VIETUS  software  should  also  work  with  Windows®  8,  Windows®  10,  Excel®  2013,  and  Excel®  2016, 
but  these  versions  have  not  been  specifically  tested.  Use  of  the  operating  system/Excel  versions  specified 
above  is  encouraged;  otherwise  the  user  is  advised  to  carefully  check  results. 

4.2  Installing  and  Starting  VIETUS 

The  VIETUS  software  is  distributed  as  a  compressed  file  (in  ‘zip’  format)  containing  four  files.  The  files 
consist  of  an  Excel  workbook  file,  two  text  files,  and  an  electronic  version  of  this  report.  These  files  may 
be  extracted  into  any  convenient  directory,  but  all  three  files  must  be  in  the  same  directory  (otherwise 
VIETUS  will  not  run).  Note  that,  when  uncompressed,  the  three  files  will  occupy  a  total  of  about  90  MB 
of  hard-disk  space. 

To  run  VIETUS,  simply  open  the  Excel  workbook  file,  as  you  would  any  other  Excel  workbook.  It  is 
recommended  that  you  manage  the  use  of  the  VIETUS  software  for  different  scenarios/projects  by 
changing  the  VIETUS  workbook  file  name  (either  by  copying  and  renaming  the  original  VIETUS  file  or 
by  performing  a  “SaveAs. . .”  operation  from  within  Excel). 

The  VIETUS  workbook  relies  on  user-defined  functions  (i.e.,  macros)  to  perform  certain  calculations,  so 
macros  must  be  enabled  for  calculations  to  work  properly.  The  user  should  select  the  “Enable  Macros” 
option  when  opening  a  VIETUS  workbook.  If  not  presented  with  an  option  to  enable  macros,  the  user 
can  try  closing  and  re-opening  the  file  or  the  user  may  need  to  alter  the  macro  security  settings  within 
Excel.  Security  options  in  recent  versions  of  Excel  are  found  under  the  name  of  “Macro  Security”  or 
“Trust  Center”  (depending  on  the  version  of  Excel).  Security  options  can  be  set  to  prompt/notify  the  user 
to  confirm  whether  macros  should  be  enabled  (or  not)  for  a  specific  workbook. 
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4.3  Description  of  the  VIETUS  Workbook 


The  VIETUS  Excel  workbook  contains  three  worksheets,  including  a  software  notice,  the  “HLC” 
worksheet,  and  the  “VIETUS”  worksheet.  The  software  notice  worksheet  provides  important 
notification/disclaimer  text  relevant  to  the  use  of  the  software.  The  “HLC”  and  “VIETUS”  worksheets 
are  described  below. 


The  “HLC”  worksheet  is  a  repository  for  contaminant-specific  information,  including  molecular  weight 
[Wieser  and  Berglund,  2009],  vapor  pressure  correlation  coefficients  [Yaws  et  al.,  2009],  solubility 
correlation  coefficients  [Mackay  et  al.,  2006;  Yaws  et  al.,  2009],  and  correlation  coefficients  for  the  gas 
diffusion  coefficient  [Yaws  et  al.,  2009].  The  “HLC”  worksheet  also  has  data  for  water  density  as  a 
function  of  temperature  [HBCP,  2014],  viscosity  of  air  as  a  function  of  temperature  [Lemmon,  2014], 
molecular  weight  of  water  [Wieser  and  Berglund,  2009],  and  the  gas  constant  [Mohr  et  al.,  2012]. 
Figure  7  shows  the  list  of  contaminants,  molecular  weights,  correlation  coefficients,  and  associated 
equations  (additional  annotations  and  data  for  water  density/air  viscosity  are  not  shown).  The  user  will 
not  generally  need  to  visit  the  HLC  worksheet  because  it  is  simply  a  data  repository. 


Correl.  Coef.  for  Calculating  Vapor  Pressure 

Correlation  Coefficients  for  Calculating  Solubility 

Correl.  Coef.  for  Calculating  Gas  Diffusion  Coef. 

Contaminant 

Abbrev. 

Contaminant 

Weight 

(g/mol) 

A 

B 

C 

E 

q 

A 

B 

C 

D 

E 

F 

A 

B 

C 

CT 

Carbon  Tetrachloride 

153.823 

7.01144 

1278.54 

232.888 

2 

9.7842  E-2 

-1.4942  E-3 

3.5854  E-5 

2.2775  E-7 

0 

0 

-4.1620  E-2 

2.5749  E-4 

5.3598  E-7 

CF 

Chloroform 

119.378 

7.11148 

1232.79 

230.213 

2 

9.9500  E-1 

-1.0531  E-2 

7.9819  E-5 

6.6431  E-7 

0 

0 

-4.9960  E-2 

2.8829  E-4 

5.9407  E-7 

DCM 

Dichloromethane 

84.933 

7.11464 

1152.41 

232.442 

2 

1.9610 

-4.4883  E-2 

8.6617  E-4 

4.9463  E-6 

0 

0 

-5.6860  E-2 

3.3534  E-4 

6.8151  E-7 

CM 

Chloromethane 

50.488 

6.99771 

870.17 

235.586 

3 

1.4019 

-6.3562  E-2 

1.7198  E-3 

-2.8262  E-5 

2.5268  E-7 

-9.3470  E-10 

-6.9250  E-2 

4.2133  E-4 

8.2225  E-7 

PCE 

Tetrachloroethene 

165.833 

7.06892 

1458.45 

226.986 

1 

2.6479  E-1 

-1.5487  E-3 

2.4477  E-6 

0 

0 

0 

-3.6840  E-2 

2.3847  E-4 

5.1123  E-7 

TCE 

Trichloroethene 

131.388 

6.87981 

1157.83 

202.58 

1 

1 .4049 

-8.2223  E-3 

1.3218  E-5 

0 

0 

0 

-4.1270  E-2 

2.6580  E-4 

5.5704  E-7 

1,1-DCE 

1,1-Dichloroethene 

96.943 

7.21678 

1181.12 

240.84 

1 

6.2741 

-3.8257  E-2 

6.0461  E-5 

0 

0 

0 

-4.5030  E-2 

2.9721  E-4 

6.0085  E-7 

cDCE 

cis-1 ,2-Dichloroethene 

96.943 

7.21953 

1290.28 

236.887 

1 

2.7735  E+1 

-1.7832  E-1 

2.9328  E-4 

0 

0 

0 

-4.6220  E-2 

3.0091  E-4 

6.1981  E-7 

tDCE 

trans-1 ,2-Dichloroethene 

96.943 

7.21356 

1244.35 

239.497 

1 

7.8039 

-4.5457  E-2 

6.9676  E-4 

0 

0 

0 

-4.6380  E-2 

3.0373  E-4 

6.2044  E-7 

VC 

Vinyl  Chloride 

62.498 

6.91423 

911.15 

239.8 

4 

-4.6670  E+2 

1 .3479  E+4 

1 .8523  E+2 

-1.3236  E-1 

0 

0 

-5.4240  E-2 

3.6658  E-4 

7.2028  E-7 

1,1,1,2-TeCA 

1 , 1 ,1 ,2-Tetrachloroethane 

167.849 

7.03897 

1467.16 

222.34 

1 

2.1790 

-1.3966  E-3 

2.9328  E-5 

0 

0 

0 

-3.6290  E-2 

2.2950  E-4 

4.9507  E-7 

1,1,2,2-TeCA 

1,1,2,2-Tetrachloroethane 

167.849 

6.91043 

1378.88 

197.086 

1 

4.8798 

-3.0937  E-3 

5.2051  E-5 

0 

0 

0 

-4.0420  E-2 

2.2899  E-4 

4.9882  E-7 

1,1 ,1-TCA 

1,1,1-Trichloroethane 

133.404 

7.00718 

1253.2 

229.624 

1 

1.0909 

-6.5278  E-3 

1.1075  E-5 

0 

0 

0 

-4.5040  E-2 

2.5770  E-4 

5.3553  E-7 

1,1,2-TCA 

1,1,2-Trichloroethane 

133.404 

7.14357 

1457.65 

228.099 

1 

2.8980 

-1.8585  E-2 

3.4896  E-5 

0 

0 

0 

-4.4400  E-2 

2.5339  E-4 

5.4055  E-7 

1,1 -DCA 

1,1-Dichloroethane 

98.959 

7.18316 

1269.43 

237.755 

1 

9.4136 

-5.7249  E-2 

9.1784  E-5 

0 

0 

0 

-4.7450  E-2 

2.8979  E-4 

5.9566  E-7 

1,2-DCA 

1,2-Dichloroethane 

98.959 

7.29525 

1407.85 

235.48 

1 

1.7915  E+1 

-1.1684  E-1 

2.0003  E-4 

0 

0 

0 

-4.8920  E-2 

2.8847  E-4 

6.0315  E-7 

CE 

Chloroethane 

64.514 

7.13047 

1097.6 

246.009 

2 

5.8420  E-1 

-1.6863  E-3 

9.3949  E-5 

-2.5316  E-6 

0 

0 

-5.5930  E-2 

3.4520  E-4 

6.8934  E-7 

x  =  (xp  / 100) 

(MWw  /  MWi)  =  xf  •  (MWw  /  MWi) 

Log(Pvap)  =  A 

B/(Tc  +  C) 

Yaws  et  al.,  2009 

1: 

xp  =  A  +  B-Tk  +  C-Tk2 

Mackay  et  al.,  2006 

Da  =  A  +  B-Tk  +  C-Tk2 

Yaws,  2003 

H  =  [  (P  /  760)  /  x  ]  ■  [  MWw  /  (pw-  Rgas'Tk) 

P vap  =  Vapor  Pressure  (mm  Hg) 

2: 

xp  =  A  +  B-Tc  +  C-Tc2  +  D-Tc3 

Mackay  et  al.,  2006 

Da  =  Gas  Diffusion  Coeff.  in  Air  (cm2/s) 

3: 

xp  =  A  +  B-Tc  +  C-Tc2  +  D-Tc3  +  E-Tc4  +  F-Tc5 

Mackay  et  al.,  2006 

MW  from  Wieser  &  Berglund,  2009  | 

Tc  =  Temp,  in  °C 

Tk  =  Temp,  in  K 

4: 

Log(xf)  =  A  + 

B/Tk  +  C-Log(Tk)  +  D-Tk 

Yaws,  2003 

xp  =  mass  fraction  x  1 00%  xf  =  mass  fractionl 


Figure  7.  View  of  the  primary  data  (molecular  weights  and  correlation  coefficients)  on  the  “HLC” 
worksheet  of  the  VIETUS  workbook 


The  “VIETUS”  worksheet  is  the  primary  user  interface,  on  which  the  user  enters  data  and  the  results  for 
the  vapor  intrusion  estimate  are  displayed.  Columns  A:F  on  the  “VIETUS”  worksheet  contain  several 
blocks  of  cells  for  entry  of  input  parameters,  a  block  of  cells  displaying  intermediate  calculation  results, 
and  a  block  of  cells  for  the  final  results,  as  depicted  in  Figure  8.  By  default,  the  “VIETUS”  worksheet  has 
space  for  up  to  three  independent  site  scenarios  (columns  D,  E,  and  F).  Input  parameters  are  segregated 
into  several  blocks  (with  shading  in  light  blue  hues),  with  each  block  containing  related  parameters. 
Moving  from  top  towards  the  bottom,  the  first  block  contains  parameters  related  to  the  nature  of  the 
source  area  and  the  vadose  zone  contaminant  transport  scenario.  The  second  block  contains  primary 
vapor  intrusion  parameters  that  must  be  specified  by  the  user.  The  third  block  contains  primary  vapor 
intrusion  parameters  that  can  either  be  directly  specified  by  the  user  or  which  may  be  automatically 
calculated  based  on  certain  assumptions  (discussed  in  Sections  0  and  2.5).  The  fourth  block  contains  the 
secondary  vapor  intrusion  parameters.  Additional  blocks  of  cells  show  the  calculated  intermediate  values 
(light  green  shading)  and  concentration  estimate  results  (tan  shading). 


30 


F  G  H 


o<  tor  UnMt-i 


&**r  input  -  lou»f  Tunny  ort  tnnlin 

iMMCMI  4 MM  t 

-  *  *  ■■■■■ 

/  ~  "T 

IF  |M|MMM«  to 

n  Mto**  to  1  • 

■to  ■MMutoV  l  ■  •  t* 

tMN-lMiMi  to  8H 

M9toPM»IT  1  MM  0  MM 

to  to*  to*  c  »  vmm*  m*»-  l»  HP  1  HU 

C  M»|  to*  GtoM wmmm  **•»’  Hi  1  8H 

mtottotoacn  to-  H 

=  1 

Source/Transport 
Input  Parameters  Block 

Vapor  Intrusion  Input 
Parameters  Block 

Block  of  Vapor  Intrusion  Input 
Parameters  that  are  either  directly 
specified  or  calculated 

Secondary  Vapor  (Needed  lf  values  ln  tlle 
Intrusion  Input  preceding  block  are 
Parameters  Block  not  directly  specified.) 

Calculated 
Parameters  Block 

Results  Block 


Notice  VIETUS  HLC 


EHC 


Figure  8.  Groupings  of  worksheet  cells  into  blocks  for  input  parameters  (blue  hues),  intermediate 

calculations  (green  hues),  and  results  (tan/orange  hues).  The  input  parameters  are  the  only 
items  that  the  user  will  modify. 


Figure  9  shows  the  details  of  the  “VIETUS”  worksheet,  including  parameter  descriptions,  variable  names, 
and  units  of  the  parameters.  For  user  convenience,  additional  information  is  presented  to  the  right  of  and 
below  the  data  input/calculational  blocks.  This  additional  information  includes  a  figure  representing  the 
generalized  conceptual  model  (to  clarify  the  meaning  of  parameters)  and  tables/notes  about  the 
permissible  ranges  of  parameter  values.  While  the  on-worksheet  information  is  useful  for  quick 
reference,  the  user  should  refer  to  this  document  for  details. 

One  additional  (optional)  input  is  present  to  the  right  (in  cell  L28);  this  input  affects  the  range  for 
moisture  content  and  is  discussed  in  Section  4.4,  along  with  information  about  how  to  use  the  tool. 
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A 

B 

C 

D 

E 

F 

21  VIETUS  (Vapor  Intrusion  Est.  Tool  for  Unsat-zone  Source) 

1  0  2  beta 

22  Descnbed  in 

PNNL-23381.  Vapor  Intrusion  Estimation  Tool  for  Unsaturated-Zone  Contaminant  Sources 

2016-Jan-06 

23 

24  User  Input- 

Source/Transport  Parameters 

25 

Scenano  Name 

— 

Case  A 

Case  B 

CaseC 

26 

Contaminant 

— 

CT 

1 ,1 ,1-TCA 

■  -t™. 

27  T 

Temperature 

rci 

19.595 

20 

28  O) 

Avg.  Moisture  Content 

[wt  %] 

4.8504 

0.808 

29  VZT 

Vadose  Zone  Thickness 

lm] 

60 

10 

30  LI 

Depth  to  Top  of  Source 

[m) 

225 

082 

31  z 

Source  Thickness 

[m] 

15 

1 

32  w  (=  1) 

Source  Width  (=  Length) 

[ml 

20 

10 

33  q 

GW  Darcy  Velocity 

[m/dayl 

0.3 

0.005 

5 

34 

Source  Strength  Input  Type 

— 

Gas  Concentration 

Mass  Discharge 

35  C„ 

Source  Gas  Concentration 

[ppmv] 

159 

36  Msrc 

Source  Mass  Discharge 

1 '  -■] 

100 

39 

Total  porosity  *  0.3.  dry  bulk  density  •  1  855  g/mL.  and 

recharge  *  0.4  crrVyr 

40  User  Input- 

Vapor  Intrusion  Parameters 

41  dx 

Longitudinal  Distance  to  Building 

[m] 

50 

15 

100 

42  dy 

Transverse  Distance  to  Building 

[m] 

10 

0 

50 

43  dz 

Depth  of  Basement/Foundation 

[m] 

4 

0.5 

2 

44  Lc 

Foundation/Wall  Thickness 

(m] 

0.1 

0.1 

0.2 

45  Qb 

Building  Ventilation  Rate 

fm’/dayl 

2000 

1800 

12000 

46 

Re-set  the  Calcs. 

47  User  Input- 

Vapor  Intrusion  Parameters*  (directly  specified  or  calculated  from  2°  parameters) 

for  Qs,  A.  etc.  | 

48  Q$ 

Gas  Flow  Into  Building  from  Soil 

[m*/day] 

8.0 

2287 

•“ 

49  Am|,i 

Inlet  Area  (Cracks  or  "Pores") 

lm*] 

0.6 

0.23 

320 

50  0C, total 

Total  Porosity  within  Crack/Fndn 

H 

0.3 

0.3 

0.3 

51  0c,g  Gas-Filled  Porosity  for  Cracks/Fndn 

H 

0.225 

0.225 

0.225 

52 

53  User  Input  -  Secondary  Vapor  Intrusion  Parameters 

54 

55 

56 

57 

58 

59 
60 
61 
62 

63 

64 

65  Calculated  Parameters/Intermediate  Values 


Vapor  Intrusion  Path  Type 

— 

Cracks 

wb 

Building  Width 

[m] 

10 

U 

Building  Length 

[m] 

20 

AP 

Building  Underpressure 

[Pa| 

kv 

Near-Building  Soil  Gas  Permeability 

(rrfl 

Ac  from  Crack  Ratio/Width/Direct: 

— 

Width 

n 

Crack  Ratio  (Ac  /  Ac): 

H 

wc 

Crack  Width 

[m] 

0.01 

Xc 

Cumulative  Lateral  Crack  Length 

[m] 

60 

Ac 

Total  Area  of  Foundation  Cracks 

[m2] 

0.6 

Porous  Foundation 
10 
20 


68  STR  Source  Thickness  Ratio*:  [-]  0.250 

70  SA  Areal  Footprint  of  Source*:  [rrf]  400 

72  RSP  Relative  Source  Position*:  [-]  1.00 

73  L 2  Distance  from  Source  Bottom  to  GW:  [m]  22  50 

74  H  Henry's  Law  Constant”:  [-]  0.890 

75  Da  Gas  Diffusion  Coef.  in  Air”:  [m*/day]  0.6885 

76  Dc  Gas  Diffusion  for  Cracks/Fndn:  [mJ/day]  0.0530 

88 

89  Results  -  Estimated  Contaminant  Concentration  at  Selected  Location  &  in  Building  (plus  Groundwater) 


0.100 
100 
0  10 
8  18 
0.550 
0.6612 
0.0509 


#DIV/0! 

0 

#DIV/0! 

0.00 

#N/A 

#N/A 

#N/A 


94 

C„ 

Final  Soil  Gas  Concentration: 

[ppbv] 

1487 

446 

#N/A 

95 

C. 

Building  Gas  Concentration: 

[ppbv] 

5.95 

0.57 

#N/A 

96 

Cw 

GW  Conc'n  Below  Bldg 

[pg/L] 

66 

252 

#N/A 

97 

98 

99 

100 
101 
102 

103  * 

104 


*  See  table  at  right  for  permissible  ranges  of 
calculated  values  (STR.  SA.  RSP  and  L2). 

*  See  the  'HLC  worksheet  for  details  of  the 
temperature-dependent  calculation  of  H  and  D*. 

*  Enter  the  value  directly,  specify  missing  1°  or  2° 
parameter  values,  or  adjust  invalid  values 

1 05  *  Basis  for  defaults  calculated  from  or  for  2“  inputs: 

1 06  Qs  for  cracks  is  calculated  per  Eqn.  IS  of 

107  PNNL-23381,  Rev  1 

1 08  Atfiitt  is  either  a  specified  crack  area  (Ac)  or  the 

1 09  surface  area  of  the  subsurface  foundation  (A,). 

110  0c,tom  is  the  same  as  for  the  vadose  zone  soil. 

111  0c.j  is  75%  Of  0e.ie«*l 

112  Xc  is  the  building  perimeter  (2U  *  2Wb). 

113  Ac  is  either  Xc  wc  or  q  Ac. 


□  Input- Primary Parameter  c°'°r  sl»dl"9  Le9e"d 
0  Input  -  Primary  Parameter  (from  manual  input  or  default  calc. ) 

0  Input  -  Secondary  Parameter 

n  Input  -  Secondary  Parameter  (from  manual  input  or  default  calc. ) 
0  Intermediate  Calculation 
0  Result  -  Intermediate/Unscaled 

□  Result -Final 

0  1°  Input  Parameter  Value  is  Not  Yet  Specified 
0  2*  Input  Parameter  Value  is  Not  Yet  Specified 

■  Parameter  is  Not  Needed/Used 

0  See  Footnote  'a'  (w  is  in  the  extended  range) 

0  Parameter  Value  is  Outside  Suggested  Range.  But  Calculations 
Will  Proceed 

0  1  *  Input  Parameter  Value  is  Outside  Permitted  Range  or  an  Invalid 
Combination  of  Parameter  Values  is  Used  (see  footnote  'e') 

■  Error  in  Intermediate  Calculation  or  Intermediate  Value  is 

Outside  Permitted  Range 

B  Error  in  Final  Result  (due  to  input  problem  or  intermediate 
calculation  error) 


Parameter 

Name 

Permissible  Range 

Key  Values 

T 

10-30 

20 

w 

1-9* 

1.5.9* 

VZT 

10-60 

10.  30.  60 

LI 

vanes” 

— 

z 

varies1 

— 

w 

10-50° 

- • 

0 

0  005  -  03 

0  005. 0  03. 0  3 

c. 

01-2000 

159 

0.1-5000 

Horn  STOMP 
»  3  morrr*  (Korea  trie 

5WWgH«Q« 


LI 

/ 

VZT 

L2 

T 

(Contaminant,  Cgs  or  Mtlc) 


00001-0002 


transverse  (relative  to  GW  low)  distance  from  source 


specified  w  V4m  ncrements 


typical  residential  bldg  -  0  1 


typcal  residential  bldg  =  1500-2400 


typical  value  =  2-14  or  <  3%  ot  Op 


equals  A.  (for  cracks)  or  foundation  surface  area 


need  enough  gas-filled  porosity  to  allow  gas  migration 


Parameter 

Name 

Permissible 

Range 

Key  Values 

STR 

0.1  -0.5 

0.1.  0.25.  0.5 

SA 

100  -  2500 

100.  400, 900.  2500 

RSP 

0.1-10 

0.1. 1. 10 

L2 

0.5-49° 

— 

H 

content  nanl-specific 

0.89 

D , 

content  nant-spec  ific 

— 

Dc 

contaminant-specific 

— 

'  The  pre-modeled  scenarios  actually  use  residual  saturation  (S,).  not 
gravimetric  moisture  content  However,  gravimetric  moisture  content  is 
requested  as  the  input  parameter  for  user  convenience  The  key  values  for  S, 
were  0.05.  0.3.  and  0.55,  which  correspond  to  moisture  content  values  of 
0.8071721,  4  843.  and  8  87889,  respectively,  for  a  subsurface  temperature  of 
20  °C.  To  keep  the  informdion  condse.  the  permissible  moisture  content 
range  is  specified  here  as  1  to  9  vtf%.  Values  at  or  above  8  87889  wt%  are 
treated  as  S,  values  of  0.55.  If  the  user  apples  the  ‘Alow  w  down  to  Sr  = 
0  05?"  option  (located  in  cell  L30.  to  the  right  of  the  table),  then  w  values 
down  to  0  8071721  wt%  (at  20  SC)  are  allowed  Otherwise,  w  values  below 
1  wt%  result  in  an  incomplete  calculation 
!  The  LI  and  L2  ranges  are  variable  (with  a  max  range  of  0.5  -  49  m)  because 
they  are  a  function  of  the  z  and  VZT  input  values  and  the  permissfcle  range 
for  RSP 

:  The  range  for  z  is  vanable  (with  a  maximum  range  of  1  -  30  m)  because  it  is  a 
function  of  the  permissible  range  for  STR  and  the  input  value  of  VZT. 

1  The  range  for  w  is  a  function  of  the  permissible  range  for  SA  and  the  source 
area  footprint  SA  key  values  are  listed  in  the  table  for  calculated  parameters 
'  The  combination  of  the  foundation/wall  thickness.  Le.  the  crack/foundation 
area,  A**,  the  gas  drffusion  coefficient  for  the  cracks/foundation,  De,  and  the 
Q,  value  must  not  cause  the  exponential  function  to  resul  in  an  error  in  the 
follcwing  expression:  expftU  QJ/fA,^  DJ], 

The  upper  limit  shown  for  Q6  is  a  soft  limit  with  a  hard  limit  applied  at  a  value 
of  80.000  m*/d  The  upper  limit  shown  for  Qs  is  a  soft  limit  with  a  hard  limit 
applied  at  a  value  of  0.06  Qb  (which  is  twice  the  observed  top  end  sub-siabto 
indoor  air  attenuation  factor  [EPA  2012]). 


HR _  _ . _  _  _ 
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Figure  9.  View  of  the  “VIETUS”  worksheet,  showing  details  of  the  input  parameters,  intermediate 

calculations,  results,  and  reference  information.  Scenario  3  (Test  C)  has  little  data,  showing 
the  coloration  of  empty  cells  and  cells  with  errors  (due  to  missing  and  out-of-range  data). 
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4.4  Using  the  Software 

The  VIETUS  software  is  implemented  as  an  Excel  workbook,  thus  the  standard  features  of  Excel  apply. 
Hence,  the  user  should  be  familiar  with  using  Excel.  Once  a  VIETUS  workbook  is  opened  (and  macros 
are  enabled,  as  discussed  in  Section  4.2),  performing  calculations  is  as  simple  as  entering  the  data 
required  to  describe  the  site  of  interest  using  the  context  of  the  generalized  conceptual  site  model  (Section 
2.1).  On  input  of  valid  data,  results  are  displayed  immediately. 

The  VIETUS  workbook  includes  two  features  to  help  maintain  the  integrity  of  the  calculations.  The 
associated  macro  code  is  locked  for  viewing  or  editing.  In  addition,  the  worksheets  are  protected  and  data 
entry  is  only  allowed  in  appropriate  data  input  cells. 

For  most  items,  data  entry  consists  of  entering  parameter  values  for  a  particular  scenario  (in  columns  D, 
E,  and/or  F).  Empty  data  input  cells  are  all  shaded  a  (darker)  blue  color,  with  the  shading  color  changing 
upon  entry  of  data  (to  either  a  light  blue/aqua,  light  yellow,  or  light  red,  as  shown  in  Figure  9).  Light  red 
shading  of  an  input  parameter  value  or  an  intermediate  calculation  value  indicates  that  invalid  data  has 
been  entered  for  one  or  more  parameters.  Light  yellow  shading  indicates  that  a  recommended  bound  has 
been  exceeded.  The  primary  cause  of  errors  is  likely  to  be  data  values  outside  the  permissible  ranges  or 
values  that  are  inconsistent  with  each  other,  such  as  a  source  thickness  that  is  inconsistent  with  the  vadose 
zone  thickness  and/or  LI  parameter.  The  tables  of  permissible  ranges  and  a  diagram  of  the  generalized 
conceptual  site  model  are  included  directly  on  the  worksheet  to  help  the  user  identify  issues  with  input 
data. 

Selection  lists  are  provided  for  the  contaminant  name,  source  strength  input  type,  depth  of 
basement/foundation,  vapor  intrusion  path  type  inputs,  and  crack  area  calculation  approach  to  ensure 
entry  of  valid  data.  A  selection  list  is  activated  by  selecting  the  input  cell  on  the  spreadsheet  then  clicking 
on  the  arrow  button  that  appears.  The  selection  of  the  source  strength  input  type  modifies  the  requested 
input  data  to  be  either  source  gas  concentration  or  source  mass  discharge,  while  graying  out  the  unused 
parameter.  Similarly,  the  selection  of  a  porous  foundation  as  the  vapor-intrusion  path  type  causes  the 
total  area  of  foundation  cracks  parameter  to  be  grayed  out.  Also,  the  selection  of  Ratio,  Width,  or  Direct 
as  the  approach  for  estimating  crack  area  will  alter  the  required  inputs. 

An  optional  input  (“Allow  ©  down  to  Sr  =  0.057”)  exists  in  cell  L28  to  indicate  whether  soil  moisture  (co) 
values  down  to  the  equivalent  residual  saturation  (, Sr )  of  0.05  should  be  allowed  (true)  or  not  (false).  This 
effectively  alters  the  permissible  range,  as  described  in  Section  2.1  (and  Figure  5).  Soil  moisture  values 
that  are  below  the  listed  lower  range  and  those  that  are  above  the  maximum  residual  saturation  are 
colored  light  lavender  to  indicate  that  special  circumstances  are  in  effect. 

The  contaminant  source/transport  input  parameters  are  relatively  straightforward,  but  the  vapor  intrusion 
input  parameters  encompass  three  blocks  of  related  inputs,  as  discussed  in  Section  4.3.  The  first  block 
consists  of  key  parameters  that  must  be  specified  by  the  user.  The  second  block  consists  of  key 
parameters  for  which  a  formula  is  provided  to  calculate  the  value  (based  on  secondary  input  parameters  or 
porous  media  properties).  The  third  block  consists  of  the  optional  secondary  input  parameters.  The  user 
has  two  options  for  setting  the  values  of  the  second  block  of  vapor  intrusion  input  parameters.  The  user 
can  enter  the  input  parameter  values  directly  (and  leave  the  secondary  vapor  intrusion  input  parameters 
unspecified).  Alternatively,  the  user  can  specify  the  secondary  vapor  intrusion  input  parameters  and 
allow  the  second  block  of  vapor  intrusion  input  parameters  to  be  calculated  automatically  based  on  the 
approaches  and  assumptions  discussed  in  Sections  0  and  2.5.  If  the  user  directly  specifies  the  value  for  a 
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vapor  intrusion  parameter  in  the  second  block  and  later  decides  to  instead  use  the  secondary  input 
parameters  to  calculate  that  parameter,  the  formulas  for  the  second  block  of  vapor  intrusion  parameters 
can  be  re-set  by  clicking  on  the  “Reset  Calcs. . button  at  row  46. 

The  final  output  consists  of  the  estimated  values  for  the  soil-gas  concentration  underneath  the  building 
and  the  contaminant  concentration  inside  the  building  due  to  vapor  intrusion.  For  reference,  the 
concentration  of  contaminant  in  the  groundwater  directly  below  the  specified  point  of  interest  is  also 
shown.  This  groundwater  concentration  arises  solely  from  the  specified  vadose  zone  source  and  the 
resulting  steady  state  contaminant  distribution.  Other  potential  groundwater  contamination  (from 
historical  and/or  non-vadose  zone  sources)  is  not  reflected  in  the  groundwater  concentration  displayed  in 
VIETUS.  The  displayed  groundwater  concentration  is  an  average  estimated  concentration  for  the  top 
meter  of  aquifer  (which  differs  from  SVEET  [Truex  et  al.,  2013],  where  groundwater  concentrations  are 
averaged  over  a  specified  well  screen  length).  Sites  with  groundwater  sources  or  with  higher  recharge 
may  have  different  groundwater  concentrations  and  may  need  a  full  site-specific  numerical  model  to 
assess  the  impacts  beyond  the  current  vadose  zone  source. 
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